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PREFACE 


In  January  1985  the  Director  of  the  Science  and  Technology  Directorate  of  the 
Strategic  Defense  Initiative  Organization  (SDIO)  asked  members  of  the  IDA  research  staff 
to  investigate  the  feasibility  of  developing  a  y-ray  laser.  The  staff  determined  what  work 
had  been  done,  who  was  currently  working  in  the  field,  and  what  work  should  be 
encouraged  or  supported.  A  workshop  was  convened  for  researchers  directly  involved 
both  in  gamma-ray  laser  work  and  in  ancillary  fields  such  as  nuclear  structure,  radiation 
propagation  in  crystals,  Mbssbauer  Effect,  and  optical  lasers.  Next,  an  in-house  study  was 
^lnd<»^taken  to  clarify  critical  issues  concerning  the  various  ptunping  schemes  proposed  at 
the  workshop  as  well  as  systems  questions  about  the  y-ny  laser  as  a  working  device. 

The  proceedings  of  the  workshop  were  published  in  the  form  of  a  report  to  the 
Iimovative  Science  and  Technology  Office  (1ST)  of  the  SDIO.  The  work  completed  in 
1985  was  presented  in  an  IDA  Paper  (Ref.  1). 

In  1986,  the  in-house  work  focused  on  extending  the  data  base,  the  nature  of 
superradiance  in  the  Y’^y  laser  context,  and  a  detailed  investigation  of  the  upconversion 
pumping  scheme.  A  discussion  of  nuclear  systemadcs,  investigations  of  electron-nuclear- 
driven  pumping,  and  lifetime  measurements  rounded  out  that  study.  The  results  of  the 
FY  1986  effort  were  also  presented  in  an  earlier  IDA  paper  (Ref.  2). 

In  1987,  the  IDA  research  staff  looked  at  the  state  of  the  art  and  assessed  the 
situation  in  y-ray  laser  work  with  focus  on  two  areas  of  research  interest  critical  to  ctxicepts 
for  developing  a  y-ray  laser.  Heating  effects  associated  with  upconversion  techniques  were 
discussed.  The  sources  of  inhomogeneous  broadening  which  destroy  the  Mdssbauer 
Effect  were  investigated  and  techniques  available  for  restoring  the  resonance  by  externa  ',  or 
internal  fields  were  considered. 

This  report  covers  the  work  done  in  FY  1988.  In  1988  the  in-house  work 
concentrated  on  (1)  establishing  a  theoretical  model  of  nuclear  superfluorescence  which 
took  into  account  specific  characteristics  of  nuclear  radiation  emission  and  transport  that  are 
important  for  superfluorescence  on  the  nuclear  scale,  (2)  examining  in  more  detail  and  for 
more  realistic  parameters  the  heating  effects  inherent  in  upconversion  concepts,  and 
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(3)  examining  the  speed  of  response  of  nuclei  to  applied  external  fields,  as  exhibited  by 
their  spectra. 

The  development  of  a  y-ray  laser  is  viewed  as  a  high-risk/high-payoff  undertaking. 
IDA'S  involvement  focuses  on  minimizing  the  risk  and  on  striving  to  redirect  the  effort 
when  proposed  schemes  are  shown  not  to  be  feasible. 
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EXECUTIVE  SUMMARY 


This  report  describes  the  1988  research  effort  by  members  of  the  IDA  staff  in  the 
field  of  Y-ray  lasers.  The  work  is  part  of  a  continuing  task  in  support  of  the  Innovative 
Science  and  Technology  Office  (1ST)  of  the  Strategic  Defense  Initiative  Organization 
(SDIO).  The  development  of  a  y-ray  laser  is  a  high-risk  science  and  technology 
undertaking.  IDA  involvement  has  focused  in  large  measure  on  minimizing  the  risk  and 
attempting  to  redirect  the  program  as  quickly  as  possible  when  proposed  schemes  prove 
infeasible.  The  report  is  presented  in  six  independent  chapters. 

Chapter  I  discusses  the  nuclear  superfluorescence  (SF)  problem.  Seven  areas  of 
concern  critical  to  the  observation  of  nuclear  SF  but  not  relevant  or  as  important  to  SF  at 
longer  wavelengths  are  discussed.  A  model  of  nuclear  SF  based  on  the  Maxwell  Bloch 
formalism  as  adapted  to  SF  by  Haake  and  Riebold  is  developed.  A  closed  form  solution  to 
the  linear  part  of  the  problem  is  obtained  and  a  direct  solution  for  the  nonlinear  part  for  the 
average  intensity  is  discussed.  The  importance  of  this  result  is  that  it  avoids  the  use  of 
Monte  Carlo  techniques  and  explicitly  exhibits  the  dependence  of  the  pulse  intensity  on  the 
strength  of  competing  processes,  (e.g.,  internal  conversion),  the  beam  attenuation  due  to 
scattering  or  absorption  processes,  (e.g.,  the  photoelectric  effect),  and  the  finite  pumping 
time  of  the  lasing  level.  The  effects  of  beam  attenuation  and  competing  transitions  on  SF 
emission  are  investigated  and  the  different  roles  attenuation  plays  in  SF  and  stimulated 
emission  are  discussed. 

Chapter  II  discusses  SF  at  wavelengths  shoner  than  the  intemuclear  distances.  A 
condition  for  SF  under  these  conditions  for  a  crystalline  solid  is  obtained  from  the 
Bonifacio-Lugiato  theory. 

Chapter  III  discusses  our  recent  investigation  of  heating  triggered  by  the 
photoelectric  effect  during  upconversion  in  a  y-ray  laser.  This  is  a  continuation  c*"  last 
year's  effort,  as  presented  in  the  1987  repon.  We  find  that  pumping  transinons  with 
energy  between  10  eV  (the  ionization  threshold)  and  10  keV  are  infeasible  because  of  this 
heating.  This  is  true  even  if  light  low  photoelectric  cross  section  materials  are  used  as  a 
heat  shield.  We  remain  skeptical  about  possibilities  for  pumping  transitions  with  energy 
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greater  than  10  keV.  Below  10  eV  no  nuclear  states  are  known.  A  proposal  to  search  for 
such  states  is  presented  in  a  subsequent  contribution. 

Chapter  IV  discusses  simple  and  inexpensive  ways  of  searching  for  nuclear  states 
located  close  to  isomeric  levels.  We  have  shown  (1987  report  and  Chapter  III  of  this 
report)  that  photoelectric  heating  eliminates  the  possibility  of  pumping  transitions  in  the 
10  eV  to  10  keV  region.  No  nuclear  states  are  known  to  lie  within  10  eV  of  isomers.  We 
propose  to  search  for  such  states  using  blackbody  radiation.  Our  calculations  indicates  that 
such  a  state,  if  connected  to  an  isomer  by  a  reasonably  strong  Ml  matrix  element,  could  be 
identified  if  it  were  about  0.1  eV  to  10  eV  above  the  isomer.  Below  0.1  eV  even  a  strong 
single-particle  Ml  transition  would  probably  have  too  small  a  cross  section  to  allow  easy 
detection.  In  the  near  future,  with  the  advent  of  tunable  coherent  sources,  the  energy 
regime  within  0.1  eV  of  an  isomer  could  also  be  searched  for  states. 

Chapter  V  describes  a  numerical  techiuque  for  solving  the  heat  diffusion  problem 
for  a  multimedia  system  when  distributed  sources  are  present.  The  technique  is  used  to 
study  the  dynamics  of  heat  removal  from  an  active  regime  by  a  heat  sink  when  the  system 
is  subjected  to  an  upconversion  laser  pulse. 

Chtq)ter  VI  calculates  the  response  of  the  nuclear  system  to  externally  varying  fields 
as  exhibited  in  the  nuclear  spectrum.  This  calculation  addresses  the  length  of  time  a  system 
has  to  be  subjected  to  line-narrowing  pulses  in  order  to  narrow  inhomogeneously 
broadened  lines  to  acceptable  widths  for  lasing. 
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I.  NUCLEAR  SUPERFLUORESCENCE 


A.  INTRODUCTION 

Siq)erfluoiescence  (and  supeiradiance)*  is  the  cooperative  spontaneous  emission 
identical  atoms  (or  other  radiators)  of  coherent  radiation  characterized  by  directed  intense 
pulses  of  duration  much  shorter  than  the  spontaneous  emission  lifetime  of  individual 
radiators  (Ref.  3). 

Generally  speaking,  the  emission  of  a  superfluorescent  (SF)  pulse  requires  the 
preparatimi  of  an  inverted  population  of  the  identical  radiatcnrs.  Experimoitally,  the  emitted 
pulses  are  characterized  by  a  SF  time  tr  and  a  delay  time  tD,  following  the  inversion  of  the 
population.  From  sin:q)le  theory.  Ref.  4,  a  relationship  can  be  obtained  between  tr  and  td> 
the  density  of  cooperating  radiators  p,  the  natural  radiative  lifetime  Xq,  the  wavelength  of 
the  emitted  light,  X,  the  number  N  of  radiatmrs,  and  the  cavity  length  /.  In  terms  of  those 
parameters^. 


— r  ’ 

fa(N), 

pX  / 

and  the  theory  predicts  a  pulse  intensity  Is(t)  given  by  the  "inverted  pendulum  expression" 


We  follow  die  convention  of  Bonifacio  and  Lugiato  (Ref.  4)  who  distinguish  between  radiation  emitted 
by  coheiently  prepared  s)rstems  and  with  initially  a  macroscopic  dipole  moment,  which  they  call 
supertadunt,  and  incoherently  prqiared  systems  which  do  not  have  a  macroscopic  dipole  moment 
initially  but  interact  through  noimal  fluorescent  decay  to  evolve  a  macroscopic  dipole  which  then 
radiates  coherently  in  a  cooperative  mode,  which  they  call  superfluorescent.  In  both  cases,  the 
intensity  is  proportional  to  N^. 

The  value  of  the  characteristic  SF  time  tR  used  here  is  that  used  by  Bonifacio  and  Lugiato  (Ref.  4)  and 

results  fipom  assuming  an  isotropic  distribution  of  the  dipole  moment  It  differs  by  a  factor  1/3  from 
the  value  derived  by  Foldm’  et  al.  (Ref.  44)  who  take  into  account  the  anisotropy  of  the  two  level 
transitioit  We  will  use  this  value  tgR  » 1/3  tr  in  further  discussions. 
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(Ref.  4),  where  t  is  measured  from  the  time  of  the  (instantaneous)  inversion. 

Qearly,  to  describe  more  realistic  experimental  conditions,  more  general  theories 
have  to  be  involved. 

It  is  expected  that  if  a  "gamma-ray  laser"  is  developed  it  most  probably  will  emit  in 
a  superradiant  (SR)orsuperfluorescent  (SF)  mode  instead  of  a  stimulated  emission  (SE) 
mode  (Ref.  S)^.  There  are  some  characteristic  features  of  superfluorescent  decay  which 
may  provide  ways  of  overcoming  specific  nuclear  problems  (Ref.  6)  which  at  present 
prevent  CW  or  pulsed  laser  (SE)  action.  The  purpose  of  this  paper  is  to  investigate  the 
possibility  of  nuclear  superfluorescence  and  describe  the  interrelationship  of  the  pertinent 
nuclear  and  solid-state  parameters  that  govern  the  experimental  realization  of  this 
phenomena.  We  also  contrast  the  nuclear  problem  with  the  atomic  and  molecular  SF 
problem. 

Superradiance  has  been  observed  in  atomic  and  molecular  systems  and  the 
phenomena  have  been  explained  theoretically  (Ref.  7).  An  excellent  discussion  of  the 
experimental  results  is  presented  by  Q.RF.  Vrehen  and  H.M.  Gibbs  (Ref.  7)  and  of  the 
present  state  of  theoretical  understanding  by  M.  Gross  and  S.  Haroche  (Ref.  8)  and 
M.F.R  Schuurmans,  Q JI.F.  Vrehen  and  D.  Polder  (Ref.  9).  Just  as  the  observation  of 
resonance  in  nuclear  systems  requires  special  conditions  that  are  nxne  stringent  than  in 
atomic  systems,  the  realization  of  superradiance  on  the  nuclear  level  must  overcome 
restrictions  that  are  often  not  important  in  atomic  and  molecular  systems.  Some  special 
consideration  has  to  be  given  to  the  short  wavelength  of  the  nuclear  radiation,  the 
crystalline  structure  and  its  effect  on  SR  and  SF,  attenuation  of  the  beam  due  to  inelastic 
scattering,  destruction  of  resonance  due  to  recoil  of  the  nucleus,  inhomogeneous 
broadening,  and  relaxation  effects.  A  strong  Mdssbauer  Effect  (recoilless  emission)  and  a 
strong  Borrmann  Effect  (enhanced  photon  transmission  in  certain  directicHis  in  a  crystal)  are 
generally  considered  indispensible  for  nuclear  superfluorescence  but  not  for  atomic  or 
molecular  superfluorescence. 

^  The  following  example  illustrates  the  reason.  An  active  cylindrical  volume,  with  a  radius  s  cm, 
and  length  »  10-2  composed  of  Fe^^  (density  »  8  x  10^2  cm-3)  will  contain  2.5  x  10^^  nuclei.  If 
the  nuclei  are  all  excited  to  an  energy  level  of  lO^ev  (lA)  with  a  natural  linewidth  of  lO'^ev, 
cofieqionding  to  a  lifetime  of  6.6  x  10-2  sec,  the  superfluorescence  time  tsf  will  be  6.9  x  10-^*  sec. 
Since  each  photon  has  the  ener^  lO^ev.  the  total  energy  emitted  will  be  2.S  x  lO^'^ev  if  the  recoilless 
fraction  is  1/10.  This  is  4  x  10-^  J;  therefore,  the  power  in  the  emitted  superfluorescence  pulse  will  be 
S.8  X  10^  W.  However,  if  the  emitted  pulse  width  were  that  of  the  natural  lifetime  the  power  would 
onlybe61W. 
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The  extensive  literature  on  superradiance  and  superfluorescence  offers  a  detailed 
physical  understanding  of  the  phenomenon  (sec,  for  example,  review  paper.  Ref.  8). 
However,  it  is  not  yet  clear  whether  any  of  the  proposed  mathematical  models  of  SR  and 
SF  can  give  reliable  quantitative  predictions  of  effects  in  nuclear  rather  than  atomic 
emission  levels. 

The  model  presented  by  R.  Bonifacio  and  L. A.  Lugiato  (B-L)  in  Ref.  4  appears  to 
rely  on  the  fewest  ad  hoc  assumptions.  It  includes  (non  relativistic)  quantum  effects,  line¬ 
broadening,  and  (non-Maikovian)  stimulation  of  the  atomic  system  by  the  spontaneously 
emitted  photons.  The  only  significant  ad  hoc  assumption  in  the  model  appears  to  be  its 
restriction  of  the  electromagnetic  field  to  a  pair  of  independent  resonant  modes.  The 
authors  partially  justify  this  assumption  by  limiting  the  geometrical  configuration  to  a 
needle-shaped  cavity  which,  if  thin  enough,  will  support  just  two  identical  waves 
propagating  in  qpposite  directions  parallel  to  the  cavity  axis  (endfiie  modes). 

Besides  the  well-known  weak  points  of  the  B-L  theory,  such  as  the  assumptions  of 
two  modes  and  the  independence  of  those  modes,  the  theory  does  not  easily  allow  for  the 
calculation  of  effects  due  to  competing  transitions  (internal  conversion,  emission  with 
recoil,  etc.)  and  transport  effects  (photoelectric  absorption).  Other  theories  have  been 
developed  which  are  based  on  the  Maxwell-Bloch  equations  and  allow  for  the  inclusion  of 
quantum  initiation  statistics  and  modeling  of  the  fluctuation  statistics.  This  is  not  of 
particular  importance  or  interest  to  the  y-ray  laser  problem  at  the  present  state  of 
sophistication  and  development  What  is  of  interest  is  that  these  theories  can  deal  in  a 
straightforward  way  with  the  phenomena  of  competing  transitions  and  photon  transport  in 
the  medium;  two  problems  of  great  sigitificance  to  tiie  y-ray  laser  feasibility  study. 

The  detailed  characteristics  of  nuclear  inversion  pumping,  radiation  emission,  and 
transport  through  media  have  an  impact  on  the  generation  of  nuclear  superradiant  pulses. 
These  phenomena  must  be  carefully  treated  in  an  analysis  of  the  feasibility  of  nuclear 
superfluorescence.  We  are  not  interested  in  the  details  of  specific  y-ray  laser  or 
superradiator  concepts  at  this  stage,  but  only  in  general  features  exhibiting  the  differences 
between  nuclear  and  atomic  systems  and  we  concentrate  on  those  that  may  not  have  been 
adequately  addressed  by  woricers  ccmsidering  atontic  systems. 

In  particular,  we  feel  that  the  following  special  issues  that  differentiate  nuclear 
superfluorescence  from  the  well-studied  atontic  supenadiance  and  superfluorescence  have 
to  be  addressed. 
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1 .  Short  Wavelength  of  the  Emitted  Radiation 

Atomic  transition  energies  are  on  the  order  of  1  eV  to  10  eV,  and  wavelengths  on  ^ 

the  order  of  10^  A  to  10^  A.  On  the  other  hand,  nuclear  transitions  of  interest  are  on  the 
order  of  10^  eV  to  10^  eV,  with  corresponding  wavelengths  on  the  order  of  about 
0.1  A  to  10  A.  The  original  Dicke  theory  of  superradiance  deals  with  systems  in  which  the 
wavelength  of  the  radiation  is  greater  than  the  diameter  of  the  volume  containing  the  ^ 

individual  radiators.  It  has  been  generalized  to  deal  with  a  less  restrictive  case:  the 
wavelength  less  than  the  diameter  of  the  volume  attaining  the  radiators  cooperating  in  the 
emission  but  greater  than  the  interradiator  distance.  In  the  case  of  nuclear  radiators,  the 
condition  has  to  be  further  generalized,  because  the  interatomic  distance  of  the  radiators 
may  be  larger  than  the  wavelength  of  the  radiation. 

The  issue  of  the  short  wavelength  of  the  emitted  radiation  as  compared  to  the 
intemuclear  spacing  (Ref.  10)  has  been  discussed  by  Trammel.  It  is  generally  ignored  by 
other  authors  dealing  with  nuclear  superradiance  (Ref.  11),  although  it  is  certainly 
controversial  as  exetr^lified  by  the  series  of  lettm  to  the  editor  following  an  article  on  the 
subject  by  H.  Lipkin  (Ref.  12).  We  discuss  the  problem  of  shorter  wavelengths  in 
Chtqrter  n  of  this  report  in  the  formalism  of  the  Bonifacio-Lugiato  theory. 

2 .  Large  Recoil  Energy  Compared  to  the  Natural  Linewidth 

A  nucleus  or  atom  emitting  a  photon  will  recoil  with  energy  Er  =  Eq^/ZMc^  where 
Eq  is  the  resonance  energy,  M  the  mass  of  the  emitter  and  c  the  speed  of  light.  A  typical 
atom  with  resonance  energy  on  the  order  of  an  eV  and  an  atomic  mass  of  50  will  recoil  with 
an  energy  of  about  10~^l  eV,  which  is  much  smaller  than  the  typical  atomic  linewidth  of 
10~^  eV.  On  the  other  hand,  the  same  atom  in  a  typical  nuclear  transition  of  10^  eV  will 
recoil  with  an  energy  of  about  10"3  cV,  which  is  orders  of  magnitude  higher  than  the 
typical  nuclear  linewidth  of  10-9  cV.  The  result  is  that  while  the  emitting  atom  is  in 
resonance  with  the  other  atoms  in  the  cavity,  the  emitting  nucleus  will  be  out  of  resonance 
with  the  other  nuclei. 

Emission  with  recoil  will  be  treated  in  our  formalism  as  a  competing  process  which 
has  the  effect  of  reducing  the  excited  state  population.  We  believe  that  in  the  present  stage 
of  the  development  of  nuclear  superradiance  this  is  an  appropriate  assumption.  Under 
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certain  well  understood  conditions  (Mdssbauer  Effect)  a  nucleus  will  emit  without  recoil, 
stnmgly  rescMiate  with  other  nuclei  and  thus  alleviate  this  problem. 

3 .  High  Multipolarity 

The  transitions  of  interest  in  nuclear  superfluorescence  are  generally  of  high 
multipolarity  (Ml,  E2,  M2,  or  higher)  unlike  in  atomic  transitions  where  electric  dipole  is 
the  prevalent  mode  of  emission.  All  the  theories  of  superradiance  known  to  the  authors  up 
to  now  have  dealt  with  electric  dipole  transitions.  Complications  with  high  multipoles 
could  arise  in  cases  where  competing  transitions  occur  (Ref.  13).  Such  considerations  will 
have  to  be  introduced  in  future  more  complete  analyses.  Problems  of  high  multipolarity 
supeifluQiescent  transitions  will  not  be  discussed  in  this  paper. 

4.  Competing  Transitions 

The  effect  of  several  transitions  competing  for  the  depopulation  of  a  level  has 
recently  been  introduced  into  theoretical  calculations  for  atomic  superfluorescence  systems. 
In  nuclear  systems  internal  conversion  is  not  only  commtm  but  strongly  competes  with  the 
radiation  transition  in  depopulation  of  low-energy  nuclear  states.^  Internal  conversion 
plays  an  especially  important  role  in  low-energy  nuclear  transitions  and  has  to  be 
considered  carefully  in  discussing  nuclear  superfluorescence.  Other  competing  electro¬ 
magnetic  transitions  are  also  comnaon  in  nuclear  systems  and  can  be  treated  similarly. 

5.  Electronic  Attenuation 

For  the  cases  of  interest,  the  linear  attenuation  coefficient  p  is  much  higher  for  the 
energy  range  involved  in  the  nuclear  transiticHi  (103  eV  to  lO^  eV)  than  in  the  energy  range 
involved  in  the  atomic  transitions  (1  eV  to  10  eV).  Multiple  passes  in  an  optical  cavity  of  a 
few  centimeters  can  be  achieved  while  nuclear  nuliation  is  limited  to  a  path  length  of  less 
than  about  0.1  mm.  The  exact  relationship  of  p  and  the  population  inversion  in  the 
generation  of  nuclear  superradiance  is  very  important  and  has  to  be  carefully  investigated. 


^  TTieintcnudcooveision  coefficient  can  lypicaUy  range  from  on  the  Older  of 

10  to  1(P  in  the  range  of  energies  of  interest  (Ref.  14). 
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6.  Finite  Pumping  Times 

The  pumping  problem  for  nuclear  SF  is  severe  and  pumping  times  on  the  order  of 
the  lifetime  of  the  lasing  level  are  usually  required.  Instantaneous  pumping  to  inversion  is 
not  expected,  thus  calculations  using  this  assumption  are  not  realistic  and  have  to  be 
abandoned.  Realistic  pumping  times  are  important  and  have  to  be  included  in  the  analysis 
of  nuclear  superradiance. 

7 .  Geometrical  Considerations 

Although  for  the  convenience  of  analysis  the  acicular  geometry  Oong  cylinder)  is 
very  desirable,  it  may  not  be  realistic  in  terms  of  the  attainable  experimental  conditions. 
The  acicular  geometry  permits  the  restriction  of  analysis  to  a  single  nxxle  emission,  and 
when  the  Iresnel  number  F  (Fad^/ZX,  where  d  is  the  cavity  diameter,  /  die  length  and  X  the 
radiation  wavelength)  is  one,  guarantees  the  optimum  diffraction  limited  condition  (Ref. 
15).  To  achieve  a  condition  for  the  active  region  determined  by  F  =  1  may  be  difficult 
experimentally.  Consideration  of  the  feasibility  of  superfluorescence  in  a  more  realistic, 
three-dimensioiud  geometry  is  desired.  Some  theoretical  work  on  three  dimensional 
geometries  has  been  presented  in  recent  publications  (Refs.  16,  17).  Further  discussions 
regarding  general  three-dimensional  geometries  will  be  left  for  subsequent  publications. 

In  this  paper  we  will  discuss  the  effect  on  nuclear  superfluorescence  of 
(1)  competing  transitions  (e.g.,  internal  conversion,  emission  with  recoil,  branching 
ratios),^  (2)  transport  effects  (electroiuc  attenuation),  and  (3)  finite  pumping  times  using 
incoherent  sources.  The  question  of  die  feasiUlity  of  supenadiance  in  a  system  of  radiates^ 
with  s^>aratkMis  longer  than  the  wavelength  of  die  emitted  radiation  is  discussed  in  Chapter 
n.  The  questions  related  to  superfluoresoence  in  a  three-dimensional  active  medium  will  be 
discussed  in  a  future  publication. 

B.  THEORY 

Superfluorescence  with  competing  transitions  and  including  transport  effects  such 
as  beam  attenuation  due  to  incoherent  inelastic  processes  can  be  modeled  using  Maxwell- 
Bloch  equations  (Ref.  8)  modified  to  include  vacuum  fluctuations.  Our  discussion  follows 


^  In  this  paper  we  treat  all  these  mechanisms  in  the  same  way,  assuming  their  impact  to  be  due  to 
their  effect  on  the  reduction  of  the  inversion  only. 
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the  Haake  and  Reibold  thewy  presented  without  derivation  in  Ref.  18.  Appendix  A 
outlines  and  discusses  a  possible  derivation  of  the  original  equations.  The  analysis  in  this 
section  which  deals  with  a  modified  version  of  those  equations  for  3  rather  than  5  energy 
levels  is  new. 

Corresponding  to  the  excitation  level  diagram  in  Fig.  1,  the  (Ref.  18)  system  of 
differential  equations  satisfied  by  the  excited  population  number  operators  Ni,  the  atomic 

polarization  operators  R^,  and  photon  field  operators  l^j  is 

fN^  =  -‘yN4  +  44  1(a) 

^  N3  = « [^1^1  ^2^32  ^2^2]-  TjNj  +  -  ^4 

^  Nj  *  +  [E32^^2  ^2^2]  -  +  Fj  N3+  ^ 

I  N,  =  +  iR;i  +  - r,N,  +  r3N3  +  i(d) 

^  R32  *  (^3  '  ^2]  ^2  ’  ^1*^21'  "2  ^  ^2*  *^2  ^2  2(a,b) 

I  =  (N3  -  N,]  4  -  4r,V  1  [r3 + 1^1  R*,  + 

^  4i  “  ^1*^2+  *  2  (^1  ^2)  2(e,f) 


-  K32  =r3'"2j  =32’‘^l‘"2r  2'‘3-‘2“^2-‘»32  2(a.b) 

1 4i  =  (N3  -  N,)  4  -  44r  i  1^3 + lii  R*, + 

^  4i  “  ^1*^2+  *  2  (^1  ^2)  2(e,f) 

In  eqs.  (1-3)  the  ±  superscripts  for  the  fields  and  polarization  operators  refer  to  the  positive 

and  negative  frequency  components  of  these  quantities.  are  stochastic  noise  operators^ 

with  zoo  means  and  second  moments  tiiat  satisfy 


^  The  stochastic  noise  q)eraiOfs  are  defined  by  their  statistks.  Polder,  Shunmuns,  and  Viehen  (Ref.  9) 
derived  the  classical  noise  equivalent  sources  from  a  quantum  mechanical  model.  They  identified  three 
distinct  components  of  the  EM  field,  namely,  (1)  vacuum  fluctuations,  Evac>  (2)  Ej^  or  dipole 
contributions  from  different  atoms  and  (3)  the  self  force  or  radiative  reaction,  Ej,  which  is  the  dipole 
contribution  from  the  local  atom.  Initially,  the  only  non  zero  component  is  Ey^c.  This  sets  up 
correlations  in  the  polarizations  of  atoms  which  produces  a  non  zero  Edd-  The  Evk  wave  vectors 
corresponding  to  the  appropriate  geometry  of  the  active  region  (cavity)  set  up  the  polarization 
appropriate  for  cooperative  emission  with  those  wave  vectors.  This  results  in  the  (^served  SF 
emission.  The  E^  field  gives  rise  to  the  natural  width  of  the  emission  line  only.  Haake  and  Reiboild 
(Ref.  18)  ose  this  theory  to  model  ndse  sources  which  trigger  SF.  Their  noise  source  is  due  to  the 
incoherent  pumping  mechanism  with  time  constant  yor  the  natural  ^xmtaneotts  decay  from  the  excited 
state  with  time  constant  r3.  In  the  present  development  we  use  the  noise  source,  as  described  in 
equations  4,  from  the  incoherent  pumping  mechanism  for  the  SF  tri^ering. 


>  =  o  , 

MJ  ^IJ  ^IJ  MJ  nj 

<  (x,t)  (x*,t’)  >  =  y5(x  -  x')  6(t  -  0  . 
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Figure  1.  Energy  level  diagram  for  the  five-state  system  modeled  by  eqs.  (1-3). 
The  two  competing  SF  transition  are  between  upper  state  3  and  lower  states 
1  and  2  and  are  Indicated  by  gi  and  g2.  respectively.  Level  4  populates 

level  3.  providing  the  inversion.  The  decay  rates  y,  r, ,  r,, 
and  Ti  are  given  by  .  ('^ji  *  ‘^32}  and  respectively.* 


*  Other  depopulation  transitions  are  generally  allowed.  We  indicate  only  those  of  interest  at  present. 


The  initial  values  of  the  operators,  referred  to  the  vacuum  state  for  the  electromagnetic  field, 
are  prescribed  by  their  vacuum  state  expectation  values: 

(x,0)I0>  =  e52(x,0)I0>  =  N3(x,0)I0>  =  0  , 


N4(x.0)I0>=1  . 


(5) 


In  addition,  boundary  conditions  implying  that  no  external  signals  impinge  on  the  system 
are  imposed: 

E*  (0,  t)  =  0  . 

The  objective  is  to  calculate  the  mean  radiation  intensity  I(t)  (at  the  right  end  of  the 
collection  of  atoms  or  nuclei),  where 

I(t)  =  <E-(/,t)E*(/,t)>  .  (6) 
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Eqs.  (1-3)  are  written  in  terms  of  the  retarded  time  so  that 
.  =  +  =  ^  where  f  and  a'  are  the 

normalized  time  and  space  coordinate.  The  real  time  and  space  comdinates  are  normalized 


t  1C ' 

to  the  superfluorescence  time  and  cavity  length,  respectively,  so  that  t'  =  ,  x'  = 

•  ^SF  * 

where  t"  and  x"  are  the  real  time  and  space  coordinates.  The  superfluorescence  time  is 

given  by 


3XV 


/  is  the  cavity  length,  X  is  the  wavelength  of  die  emitted  photon,  x  is  the  natural  radiative 
lifetime,  and  p  is  the  inversion  density  of  cooperating  nuclei.  Of  the  parameters  indicated 
in  Fig.  1,  tij  give  the  radiative  spontaneous  emission  lifetime  of  the  state  i  due  to  the 
transition  i  j,  with  corresponding  emission  rates  Fij  and  Fs,  F2,  and  F 1  give  the  total 
single  resonator  spontaneous  incoherent  emission  rates  from  states  3,  2,  and  1, 
respectively.  The  coupling  constants  for  transitions  3-4  2  and  3  -4  1  are  given  by  g2  and 
gl,  respectively. 

Calculations  of  superradiant  phenomena,  described  by  the  general  eqs.  (1-3)  and 
the  level  structure  of  Fig.  1,  using  Monte  Carlo  techniques,  will  be  described  in  a 
subsequent  publication.  Presently  we  are  interested  in  studying  (1)  die  effect  of  counting 
transitions  such  as  the  incoherent  decay  process  described  by  the  decay  rate  F3,  which 
removes  excited  state  nuclei  from  the  system,  (2)  the  effect  of  transport  phenomena  as 
included  in  the  model  through  the  linear  attenuation  coefficient  p  in  eqs.  (3a,  3b),  and  (3) 
the  effect  on  superfluorescence  of  die  pumping  rate  y. 

Equations  (1-3)  can  be  solved  for  some  special  cases  to  get  insight  into  the  effect  of 
these  parameters  on  the  phenomena  of  supenradiance.  With  this  intent  we  simplify  the 
problem  by  craisideiing  levels  4, 3,  and  1  only  and  drop  the  noise  source  The  relevant 
set  of  equations  can  be  separated  into  two  groups.  The  first  group  describes  the  formation 
of  superradiance  and  consists  of 
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The  second  group  of  equations  describes  the  population  dynamics  of  the  three  states  and  is 
used  to  calculate  the  inversion  function  f(x,t)  =  N3  -  Ni  , 


-3^=--yN,  8(a) 

^=-(E>;+^.R„)-r3N,-eiN,  8(b) 

3N,  (  ^  ^  ^ 

-^=+(E;R;,+^iR3,)+r3N3  .  8(c) 

Initially,  only  level  4  is  populated  so  that 

N4(0)  =  No 


N3(0)  =  Nj(0)  =  0  , 


and  since 


+  N3  +  Nj  =  No  for  all  time,  (9) 

we  get 

f(x.t)  =  2N3(t)+N4{t).No  .  (10) 

Equations  (7a,  7b),  after  dropping  the  subscripts,  can  be  combined  to  give 

Assurrung  that  f(x,t)  is  a  function  of  time  alone  and  taking  the  Laplace  transform  of  (1 1) 
with  respect  to  x  leads  to 


where 


4 11  +  s)  +  s)  -  g  f  Je  =  ^(0,  t)  +  Le  (0,  t)  +  g  C 


E(s)  =  J  e-®*  E(x)  dx 


.  (12) 


I-IO 


Taking  the  appropriate  initial  condition  that 

E(x.O)  =  0  , 

which  implies 

E(s,0)  =  0  . 

the  solution  of  eq.  (12)  is 


I 

E(s,t)  =  g  j  |(s.Oe 


y  (t-0 


V 


df  . 


(13) 


where 


t 

k(t.O  =  gJf(ndt"  . 


(14) 


Since  is  the  Laplace  transform  of  the  Bessel  function  Iq  (z/E) .  is  the 


s+yH 


Laplace  transform  of  e  ^  (2yEt  ] .  Therefore, 


E^(x.t)  =  g  J  J  lo  (2Vl^(t,t’)  (X  -  X’)  ]e  2  ■  2  ‘  ^  (X.  j.)  .  (15) 

O  O 


The  experimentally  measured  intension  of  the  superfluorescent  emission  is  given  as  the 
average  of  the  product  of  E+  and  E-  at  the  right  boundary,  /  =  1.  From  eqs.  (15)  and  (4) 
this  is  given  by 

I(t)  =  <E-(l,t)E^(l,t)> 

,  till 

6  0  0  0 


“  N  JJJJ  Io[2Vk(t.O(l  -x’)] 
o  o  o  o 

e'TUd  -X*)-  y  <‘-‘^1  [2Vk(t,t") (1  - x") ) 


e  2  ^  6(x’  -  X")  6(t’  - 1")  dx'  dt’  dx"  dt" 
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I 


=  ^  0^“"  j  j  4(2Vk(U')<l-x’)>'“'*'^-’*'  dx’  dl' 

o  o 

j  ti^MuO  dx’df  (16) 

o  o 


In  general,  including  the  nonlinear  terms  +  E*  R',  the  expression  for  k(t,t')  according 
to  eqs.  (10)  and  (14)  is  given  by 


k(t,f)  =  2  J  NjCO  dt"  +  N„  45 — -  N„(t  -  O  . 

t’ 


because  of 


N4  =  N„c‘^‘  . 


(17) 


(18) 


In  the  linear  region  of  validity  of  equations  (7)  and  (8),  ignoring  the  E+R+  +  E-R-  terms, 
the  population  density  eq.  (8)  lead  to 

t 

Nj  =  rJ  N,(0dl' 


=  2^  1) 


r-Y 


From  eqs.  (17)  and  (19)  we  get 


k(t,t’)=J  [2^3(0 +  N,(f')-N,]dt" 


e-n-.e-r.] 

(  .  \ 

[  Y  r  J  ^0 

1  7  J 

(19) 


-No(t-t*)  .  (20) 
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We  next  consider  the  modification  of  eq.  (19)  for  the  nonlinear  regime.  Defining 


S=-i<E-E%  . 


it  follows  fitxn  eq.  (7b)  that 


S=-g(R‘'E^  +  R-El+ ^E-E''  ,  (22) 

where  the  fact  that  the  and  the  E^  commute  has  been  used.  From  eqs.  (8b),  (19),  and 
(22)  it  follows  that 

N3  =  I^(e-^Le-r‘)  +  y  (S  -  ji  E-  e"^  )  dt '  . 


Now,  defining  an  average  N3  by 


1 

^3  =  J  <N3>  dx  , 


eq.  (23)  can  be  replaced  by 


whoe 


1^=  Je"’^^*"‘^J<E"E^>dxdt'  . 

O  O 

For  4  =  o  the  last  term  in  eq.  (24)  disappears  and  the  solution  to  the  problem  can  be  written 
as  an  iterative  procedure  in  terms  of  I(t),  the  intensity  at  right  side  of  cavity.  For  finite 
values  of  p  the  third  term  describes  a  loss  due  to  absorption  in  the  medium. 

The  third  term  in  eq.  (24)  will  be  small  for  small  values  of  p.  For  larger  values,  an 
approximation  in  which  the  integral  of  <E“E'*>  is  taken  to  be  proportional  to  the  radiation 
intensity  I  can  be  used.  Then  (24)  can  be  replaced  b 

^3  ~  ■  "g  f  ^  ^  • 

o  os-v 
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and  /  is  an  appropriate  length  lying  between  0  and  1.  It  has  to  be  determined  from  test 
runs  on  convergence  (/  =0.1  and  0.025  were  found  convenient  for  our  runs). 

Using  eq.  (25)  instead  of  eq.  (19)  in  the  calculation  of  k  (t ,  1 0  and  iterating  with 
(16)  gives  us  the  solution  to  the  nonlinear  problem.  Appendix  B  discusses  the  numerical 
approach  used  in  implementing  the  solution. 

C.  DISCUSSION  OF  RESULTS 

With  reference  to  eq.  (16),  we  can  analyze  the  effects  of 

(1)  Transitions  competing  with  the  SF  transition.  In  this  paper  we  have 
concentrated  on  the  first  two  effects.  The  modification  of  the  SF  pulse  due  to 
finite  pumping  time  will  be  treated  in  a  subsequent  publication.  Competing 
transitions  such  as  emission  with  recoil  characterized  by  the  recoilless  fraction/ 
and  electron  ejection  from  the  nucleus  as  characterized  by  the  internal 
conversion  coefficient  a,  as  well  as  other  processes  characterized  by  the 
branching  ratio  p,  although  quite  different  physically,  are  treated  similarly 
here.  It  is  assumed  that  their  sole  effect  is  to  reduce  the  inverted  population, 

(2)  Electronic  absorption,  as  characterized  by  the  linear  extinction  coefficient  (i, 
and 

(3)  The  fiiute  pumping  time  to  obtain  the  required  inversion. 

There  are  two  first-order  effects  expected  from  a  finite  linear  extinction  n.  One 
effect  is  clearly  expected  from  eq.  (16),  because  of  the  exponential  factor  e*2tix’  which 
attenuates  the  beam  as  it  traverses  the  medium.  This  would  be  exhibited  in  the  experimental 
results  as  a  decrease  in  the  intensity  of  the  emitted  superfluorescent  pulse.  The  other  effect 
is  a  result  of  the  suppressed  development  of  the  electric  field  as  expressed  by  eq.  7.  For 
increasing  values  of  |i,,  one  expects  that  a  longer  time  would  be  required  for  the  generation 
of  a  large  enough  E  to  initiate  the  SF  pulse.  This  would  be  observed  experimentally  as  an 
increase  in  the  time  delay,  to- 

Equation  (16)  was  used  to  calculate  the  pulse  shape  as  a  function  of  the  decay  rate  F 
and  the  linear  attenuation  coefficient  ^  for  ^  =  0.  Figure  2  shows  the  variation  in  the 
lineshape  as  F  is  varied.  Notice  the  expected  decrease  in  the  peak  intensity  and  the  shift  in 
the  peak  position  as  F  increases.  Figure  3  shows  the  variation  in  the  pulse  shape  as  a 


function  of  H-  Notice  again  the  decrease  in  intensity,  increase  in  the  delay  time,  and 
broadening  of  the  pulse  as  |i  increases. 


Figur*  2.  Superfluorescent  pulse  calculated  with  y  s  5  for  p  s  0  and 

different  values  of  r.  Curves  are  all  normalized  to  the  curve  label  (a).  The  r 
values  and  the  relative  Intensities  Im  (parentheses)  are: 

(a)  r  s  0.01  (1.0) 

(b)  r  s  0.2  (0.61) 

(c)  r  s  0.5  (0.30) 

(d)  r  8  1.0  (0.19) 


It  is  instructive  to  compare  the  effect  of  p.  on  lasing  and  superfluorescence.  These 
are  different  phenomena  that  have  been  unified  in  a  theoretical  description  (Ref.  19)  but  can 
be  differentiated  experimentally  as  observed  by  Okada  et  al.  (Ref.  20).  In  their 
experiments,  the  onset  of  superfluorescence  was  preceded  by  amplified  spontaneous 
emission  pulses  as  the  inversion  density  increased. 


0  10  20  30  40  50  60  70 

TIME  (in  units  of  Toe) 

4-26-90-3  • 

Figura  3.  Superfluorescent  pulses  calculated  using  the  program  described  In 

Appendix  B  for  7  »  so  0.001  -;V  .  B  1.0  and  different  values  of  p  as 

Indicated  next  to  the  curves  In  the  figure.  The  curves  are  ail  normalized  to  1. 

The  relative  Intensity  values  are  0.16,  0.11,  0.037,  0.0074  and  0.00073  for  • 

Increasing  values  of  p. 


Out  study  shows  that  the  effect  of  beam  attenuation  (characterized  by  the  linear 
attenuation  coefficient  p)  on  the  development  of  lasing  is  different  than  its  effect  on  the 
development  of  a  superfluorescent  pulse.  Hie  stimulation  cross  section  as,  inversion  An* 
[=  f  (x,  t)  =  N3  -  Ni  as  in  eq.  (7a)]  and  p  are  coupled  in  a  linear  relation  forming  the  gain 
coefficient  k  which  characterizes  the  lasing  phenomena.  The  Schawlow-Townes  condition 
for  lasing  demands  that  the  gain  coefficient  k  «  CsAn*  -  p  >  0  in  order  for  lasing  to  take 
place.  The  gain  in  a  system  exhibiting  amplified  spontaneous  emission  (ASE)  is  given  by 

G  =  5_J_ 
kJ 

For  k  >  0  gain  is  observed,  but  for  k  <  0  only  attenuation  of  the  beam  traversing  the 
medium  is  possible.  In  SF  inversion  An*  and  attenuation  p  interact  through  the  nonlinear 
relationship  given  in  eqs.  (1-3)  and  thus  their  awnbined  effect  on  SF  is  more  complex. 

For  SF,  inversion  (An*  >  0),  is  required  to  promote  the  development  of  the 
polarization  from  the  quantum  flucmation,  eq.  (7),  but  the  effect  of  p  on  the  development 
of  E  is  only  to  delay  the  SF  pulse  and  modify  its  shape  but  not  prevent  it  if  other  effects 
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such  as  competing  processes  do  not  destroy  the  inversion.  This  can  be  seen  from  the 
solution  of  eq.  (7b)  which  can  be  written  as 

_±  -J:  -UX  f 

E^j(x)  =  ^j(0)  e  *“  +  J  e  2  gjR3,  (y)  dy  . 

O 

For  large  p  the  first  term  is  unimportant  and  only  regions  of  the  integral  for  y  «  x 
contribute.  However,  there  is  a  gradual  decrease  and  not  a  sharp  cutoff  or  threshold. 
Further  quantitative  investigation  of  this  phenomenon  with  the  complete  theory  is 
necessary. 

The  pulse  shape  as  a  function  of  initial  inversion  and  the  linear  attenuation 
coefficient  p  for  amplified  spontaneous  emission  and  superfluorescence  are  compared  in 
Appendix  C.  These  results  were  obtained  assuming  instantaneous  inversion.  In  general, 
the  ASE  lineshape  peaks  at  the  natural  lifetime  and  only  the  relative  intensity  decreases  as 
the  initial  inversion  An*  or  p  decrease.  For  SF,  in  addition  to  these  effects,  the  lines 
broaden  and  shift  to  later  times  as  An*  and  p  decrease.  For  high  inversion  and  low 
attenuation,  the  SF  peak  is  much  narrower  than  die  natural  lifetime  x  and  appears  much 
earlier  than  the  ASE  pulses. 

The  effect  of  the  linear  attenuation  coefficient  (on  the  SF  pulse  shape)  has  been 
considered  as  arising  from  the  limitadon  on  the  effective  cavity  length  L,  and  thus  on  the 
maximum  number  of  cooperating  emitters  N  (Ref.  21).  This  effect  would  be  observed  as 
an  increase  in  the  delay  time  xd,  an  increase  in  the  width  of  the  pulse  xr  and  a  decrease  in 
the  intensity  of  the  pulse.  The  inverted  pendulum  expression  for  the  SF  pulse  could  then 
be  modified  to  take  this  into  account  by  assuming  that  N,  xr,  and  xd  are  functions  of  p, 
and  the  pulse  is  given  by 


The  effect  of  p  in  eq.  (26)  is  assumed  to  be  in  decreasing  the  active  volume  by  a  factor  of 
0 

,  where  6  is  an  od  hoc  parameter. 

It  is  interesting  to  compare  the  results  from  the  reduced  efrective  volume 
assumption  represented  by  eq.  (26)  with  the  results  obtained  with  the  modifred  Maxwell- 
Bloch  equations,  including  random  sources  [eqs.  (7a,  7b  and  8a,  8b,  8c)].  We  assume  that 
the  modified  Maxwell-Bloch  equations  model  SF  from  an  incoherently  pumped  two-state 
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system  correctly  while  eq.  (26)  provides  a  simpler  intuitive  but  ad  hoc  model.  The  SF 
pulses  obtained  from  eqs.  (7,  8)  for  five  values  of  )i.  arc  plotted  in  Fig.  3  with  all  peaks 
normalized  to  1.0  independently.  The  widths,  intensities,  and  delay  times  of  these  pulses 
are  given  in  Table  1. 


Table  1.  Pulse  Parameters  Derived  from  Figure  3 


p 

''D 

1 

W 

A 

0 

6 

0.16 

5.0 

0.80 

2 

8.5 

0.11 

6.5 

0.72 

10 

18 

0.037 

13.5 

0.50 

20 

29 

0.0074 

23.0 

0.17 

40 

55 

0.00073 

46.0 

0.034 

P  - 

attenuation  coefficient 

to  - 

delay  time 

1 

intensity 

W  « 

width 

A  • 

area  under  peak  (1  x  w) 

The  delay  time  is  increased  because  the  effective  length  L  of  the  active  region  is 

decreased  to  L'  »  ^  for  p. »  1.  In  Figure  4  we  compare  the  delay  times  Xp  calculated 
from  equations  7  and  8  (pulses  shown  in  Fig.  3)  with  those  computed  by  assuming  a 


,  .  .  .  -  K(u  + 1)  , 

reduction  m  cooperation  number  nom  an  increase  in  p  given  by  x^  =  — - in 


N„ 


p+1 


The  horizontal  arrows  in  Fig.  4  indicate  the  reduction  in  p  required  for  Xjj  =  x^  and 


values  of  0  give  the  actual  reduction  factors.  In  Figure  5  we  plot  and  as  a  function  of 

A  0 

p  to  compare  the  total  emission  calculated  from  the  equations  7  and  8  with  the  expected 


emission  if  the  total  volume  of  cooperating  nuclei  were  decreased  by  ^  ^  j.  A  reduction 
parameter  of  0  =  3.1-6.7  for  p  =  5  to  20  represents  consideration  of  an  effective  average 
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length  coiresponding  to  the  notion  that  not  all  the  photons  travel  the  total  length  of  the 
cavityJ  These  findings  compare  qualitatively  with  the  results  of  MacGillivray  and  Feld 
(Ref.  22)  who  calculated  the  effect  of  beam  attenuation  in  order  to  explain  their 
experimental  results  on  HF  gas  (see  their  Fig.  5).  Their  calculations  show  an  increase  in 
the  delay  time  and  a  decrease  in  the  beam  intoisity  as  the  attenuation  increases.  Specifically 
for  L  =  1,  =  2.5  cm-l,  and  5.0  cm-l  their  results  show  about  40  percent  and  60  percent 
increases,  respectively,  in  the  delay  times  fit)m  the  p  =  0  delay  time.  Our  results  (Fig.  5) 
show  larger  delay  time  increases  (about  twice  their  amount)  but  an  exact  quantitative 
compaiistm  was  not  expected  because  other  parameters  in  the  two  calculations  arc  different 
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-iii  -4/ 

The  reduction  parameter  of  0  =  3.1*6.7  is  obtained  from  the  relationship  e  ®  ^  ^  | 

andti-Sand20.  ^ 


0  10  20  30  40  50 

44640-1 

Figure  4.  The  delay  time  xg  (solid  line)  taken  from  Fig.  3  for  different  values 

of  \L  and  the  delay  time  ^  (dashed  line)  calculated  assuming  \i  reduces  the 
cooperative  volume  plotted  as  a  function  of  p.  The  values  of  e  represent 
reduction  factors  for  different  p  values  so  that  ^  >  v  The  circles  show  points 
calculated,  the  connecting  lines  are  drawn  to  help  the  eye  follow  trends. 
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FIgura  5.  Plot  4-  and  7  at  a  function  of  \i  for  tha  diffarant  pulaaa  ahown  In 

Fig.  3.  Again  tha  diffarant  vaiuaa  of  e  rapraaant  raductlona  factora  on  p  that 
daacriba  affactivanaaa  of  p  in  raducing  cooparation  voluma.  CIrclaa  Indicata 
caleulatad  pointa,  tha  connacting  llnaa  wara  drawn  to  halp  indicata  tranda. 
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D.  SUGGESTED  EXPERIMENTAL  INVESTIGATIONS 


Our  studies  on  pumping  of  isomeric  levels  (Ref.  23)  and  the  nuclear  SF  calculations 
discussed  here  indicate  that  experiments  examining  the  interaction  of  the  linear  extinction 
coefficient  p,  the  competing  phenomena  characterized  by  Fs,  and  finite  pumping  times 
characterized  by  Y  should  be  undertaken.  We  propose  experiments  performed  in  the  optical 
or  longer  wavelength  regime  where  SF  has  already  been  observed,  be  undertaken  to  check 
the  theory  developed  here.  Systems  composed  of  CH3F,  HF,  Na,  or  T1  could  be  used 
for  this  purpose  since  SF  in  them  has  already  been  observed.  Introduction  of  varying 
concentrations  of  inert  gases  could  be  used,  for  example,  to  study  the  dependence  of  SF  on 
p.  in  a  gaseous  medium.  Variation  in  the  length  and  intensity  of  pumping  laser  sources 
could  be  used  to  check  the  effect  of  the  pumping  rate  y,  and  the  effect  of  competing 
transitions  could  be  studied  by  controlling  the  inversion  of  different  transitions  through 
selective  population  of  the  lower  levels  (2  and  1  in  Fig.  1),  or  through  manipulation  of 
relaxation  rates  through  temperature  or  external  fields.^ 

The  second  set  of  experiments  we  propose  deal  with  nuclear  transitions  aiul  the 
observations  of  modifications  of  the  nuclear  decay  rates.  These  experiments  can  be  done 
with  the  isomers  identified  in  Ref  22.  Pumping  by  thermal  neutrons  to  produce  inversion 
in  a  symmetric  sample  would  be  followed  by  observation  of  the  total  count  rate  from  the 
isomeric  level  using  a  4ic  detector  to  catch  most  of  the  radiation.  Altiiough  non-coq)erative 
eixiission  would  probably  predominate,  and  only  properly  oriented  and  prepared  <maii 
regions  of  the  sample  may  produce  SF  pulses,  these  may  be  enough  to  be  observed 
through  decay  curve  modifications.  We  expect  these  experiments  with  arbitrary  sample 
geometries  would  be  much  easier  than  single- whisker  experiments  usually  described  in  the 
literature.  Further,  more  detailed  discussions  of  these  two  types  of  experiments  will  be 
presented  elsewhere  in  future  publications. 

E.  CONCLUSIONS 

We  have  identified  several  features  of  nuclear  transitions  and  radiation  transfer 
through  materials  which  have  to  be  incotpcxated  into  an  SF  theory  in  order  to  determine  the 
feasibility  of  nuclear  SF.  The  Haake-Reibold  model  of  SF  was  generalized  to  treat 
transport  phenomena  and  used  to  study  the  effects  of  competing  transitions,  electronic 


*  Further  diacussions  of  the  effects  will  be  left  to  future  publications. 
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attenuation  of  the  photon  beam,  and  finite  pumping  time  on  the  development  of  SF  pulses. 
It  was  determined  that  the  electronic  attenuation  slows  down  the  development  of  the  SF 
ptilse  (increases  Td)  and  lowers  the  emitted  intensity  but  does  not  produce  a  hard  threshold 
on  SF  as  it  does  on  ASE,  as  characterized  by  the  Schawlow-Townes  gain  requirement 
The  competing  transititais  like  (1)  internal  conversion,  (2)  branching  transitions  producing 
the  emission  of  photons,  and  (3)  y  emission  witii  recoil  (non-M5ssbauer  transitions)  reduce 
the  SF  intensity  but  do  not  affect  a  strong  shift  in  the  pulse  emission  time  until  the  rate  of 
these  processes  approaches  the  SF  delay  time.  All  of  these  predictions,  as  well  as  the  effect 
of  the  finite  punning  time  and  geometiy  should  be  investigated  further. 

Finally,  as  a  way  of  investigating  the  effects  of  electronic  attenuation,  competing 
transitions,  and  finite  pumping  time  on  SF,  we  propose  that  experiments  at  longer 
wavelength  be  performed  using  systems  and  conditions  under  which  SF  has  been 
observed.  Such  experiments  could  be  much  more  easily  performed  than  nuclear 
experiments  and  would  be  valuable  in  verifying  the  themetical  results  presented  here. 
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n.  SHORT  WAVELENGTH  SUPERFLUORESCENCE 


If  wavelengths  associated  with  nuclear  transitions  are  smaller  than  interatomic 
spacings,  efficient  radiation  can  only  occur  when  conditions  for  Borrmann  mode 
propagation  are  satisfied.  Since  superfluorescence  is  also  necessary,  the  effect  of 
wavelengths  shorter  than  interatomic  spacings  on  superfluorescence  is  also  of  interest.  In 
particular,  it  is  of  some  interest  to  determine  the  extent  to  which  conditions  for 
superfluorescence  at  short  wavelengths  may  be  in  conflict  with  those  for  Borrmann  mode 
propagation. 

As  observed  in  the  introduction  to  Chapter  I,  in  nuclear  superfluorescence 
wavelengths  associated  with  transitions  of  interest  will,  in  fact,  be  smaller  than  the  spacing 
of  the  contributing  nuclei.  The  behavior  of  a  similarly  spaced  classical  periodic  array  of 
coherent  sources  suggests  that  for  a  potentially  superfluorescent  crystal,  regarded  as  an 
array  of  identical  nuclei,  a  first-order  effea  of  such  a  short  wavelength  would  be  to  limit 
collective  radiation  to  certain  directions.  In  fact,  the  relatively  simple  model  of  atomic 
super-fluorescence  due  to  Bonifacio  and  Lugiato  (Refs.  4  and  24)  has  sufficient  scope  to 
ccxifiim  the  existence  of  such  an  effect 

Hie  Bonifacio-Lugiato  superfluorescent^  model  for  two-state  atoms  in  a  rectangular 
crystal  lattice,  as  elucidated  in  Refs.  4  and  24,  considers  only  the  case  in  which  the  emitted 
radiation  wavelength  is  larger  than  the  atomic  spacing.  However,  this  limitation  is  not 
inherent  in  the  model,  despite  its  simplistic  treatment  of  propagation  and  other  possible 
geometrical  effects. 

Reference  25  treats  the  case  of  nuclear  transitions,  for  which  the  wavelength  is 
smaller  than  the  atomic  spacing,  but  the  paper's  discussion  of  superfluorescence  is  less 
informative  than  that  of  Refs.  4  and  24.  This  is  partly  because  the  Bonifacio-Lugiato 
model  uses  a  representation  for  the  coUective  excitation  operator  that  relates  it  to  the  lattice 
structure,  whereas  Ref.  25  uses  a  representation  that  relates  the  operate  to  the  phase  of  the 
emitted  photons  at  the  lattice  sites. 

In  particular,  when  the  wavelength  of  the  emitted  radiation  is  less  than  the  atomic 
spacing,  the  B&L  model  implies  that  superfluorescence  can  only  occur  when  the 
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wavelength  and  the  crystal  lattice  structure  are  related  in  a  certain  way.  In  fact,  the  required 
relation  is  also  sufficient  for  coupling  of  the  emitted  radiation  to  a  Boirmann  mode. 

Reference  26  gives  an  analysis  of  superradiance  that,  like  the  B&L  model,  relates 
the  collective  excitation  to  the  lattice  structure  and  should  therefore  be  capable  of  yielding 
the  short  wavelength  superfluorescence  conditions.  Because  its  treatment  includes  more 
geometrical  detail  than  the  treatment  of  Refs.  4  and  24,  for  the  long  wavelength  case 
Ref.  26  is,  in  fact,  able  to  relate  the  spatial  pattern  of  the  radiation  emitted  by  the  collective 
atomic  array  to  the  characteristic  shape  of  the  active  volume.  The  analysis  of  Refs.  4  and 
24,  however,  is  restricted  to  the  case  in  which  the  active  volume  is  a  thin  needle. 

On  the  other  hand.  Ref.  26  shows  that  the  needle  shape  provides  the  maximum 
gain,  albeit  in  a  single  direction.  Thus,  to  the  extent  that  it  is  possible  to  implement,  the 
needle  is  presumably  the  preferred  shape  for  a  superradiant  structure.  Therefore,  it  seems 
useful  to  consider  the  consequences  of  the  B&L  model  for  the  short  wavelength  radiation 
that  characterizes  the  nuclear  case. 

The  interaction  Hamiltoiuan  originally  introduced  by  B&L  has  the  form 
N 

Hj  =  ^  5]  8k  {  -  ©j)  t  -  ik-xp  -  H.C. }  ,  (27) 

vV  j,i  k 

where  k  is  die  emitted  photon  wave  vecmr  for  a  particular  mode,  xj  is  the  position  vector  of 
the  jdi  lattice  site,  a^  is  the  corresponding  photon  creation  operator,  cuo  is  a  mean  reference 

frequency,  and  rr  is  the  de-excitation  q)etator  fw  the  atom  at  the  j*  lattice  site.  B&L  also 

define  the  reciprocal  lattice  vectors  a,  each  component  of  which  has  the  form  (27c/L)n,  n  = 
0, 1, ...,  N  -  1,  where  L  is  the  corresponding  lattice  dimension.  They  define  the  collective 
excitation  and  de-excitation  operators  R^a)  as  functions  of  a  by 

N 

(a)  =  X  ^  ”P  .  (28) 

j=i 


from  which  it  can  be  shown  that 


On  substituting  eq.  (29)  into  eq.  (27)  it  follows  that 


R"  (a)  f*  (k  -  a ,  t)  -  H.C.]  •  (30) 

VV  a  K 


I 


P 


where 


f(Tl,t)  =  i  J^exp(iiTXj)exp[i(a)j~o>o)t]  (31) 

j=i 

The  time-dependent  factors  in  eq.  (31)  are  separated  out  on  the  assumption  that  they  can  be 
replaced  by  an  average  over  all  atoms  and  that,  because  the  frequencies  coj  associated  with 
the  individual  atoms  are  uncorrelated,  they  can  be  removed  from  f  (t]  ,  t)  and  subsumed  in 
the  coupling  constant  gk,  which  then  becomes  time  dependent.  This  time  dependence  of  gk 
is  a  way  of  introducing  inhomogeneous  broadening. 

At  this  point,  B&L  introduce  the  assumption  that  the  atomic  spacing,  given  by 

d  =  UN  , 

is  much  smaller  than  the  photon  wavelength,  which  implies  that  the  function  frT])  (wherein 
the  time  dependence  is  no  longer  indicated  for  the  reason  just  noted)  is  large  when  q 
vanishes  and  is  otherwise  small.  However,  if  this  assumption  is  not  made  then  frq)  will  be 
large  whenever  q  •  xj  =  27n)j,  "Oj  =  0,  ±  1, ...  and  otherwise  small. 

The  inteiacdon  that  results  in  photon  emission  when  cooperative  de-excitatiem  of  the 
atoms  takes  place  can  only  occur  when  f(k  —  a)  in  eq.  (30)  is  not  negligible  for  some  value 
of  k  and  a,  i.e.,  when 


(k  -  a)  •  X.  =  litM.  ,  (32a) 

where  Vj  is  in  the  set;  0,  ±  1, ...,  which  implies  that 

d/X-n/N  =  \),n  =  0, 1,...,N-1  ,  (32b) 

where  v  is  in  the  set:  0,  ±  1, ....  In  eq.  (32)  it  is  tacitly  assumed  that  the  propagation 
direction  is  along  a  narrow,  needle-shaped  active  volume,  in  which  only  a  single  plane 
wave  photon  mode  is  supported.  More  generally,  X  should  be  intopreted  as  Xj  the  photon 
wavelength  divided  by  the  component  of  the  unit  vector  in  the  propagation  direction  along 
the  lattice  direction  defined  by  the  xj. 
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It  follows  from  eq.  (32b)  that  the  possible  values  of  m  must  be  non-negative.  Thus, 

=  lV(N\)j  +  np,  0<nj<N,  0  .  (33) 

On  setting  dj  =  0, 1, ...  in  succession  it  follows  from  eq.  (33)  that  for  Oj  =  0,Xj>d, 
for  Vj  =  l,d>Xj>d/2,  for  Vj  =  2,d/2>Xj>d/3,  etc.  That  is,  for  any  given  photon  wavelength 
one  and  only  one  value  of  exists  that  will  satisfy  the  required  condition. 

Referrace  25  gives  as  the  condition  for  a  Borrmann  mode  the  relation 

|k-al  =  lkl  (34) 

for  some  reciprocal  lattice  vector  a.  To  satisfy  eq.  (34)  it  is  sufficient  that  k  =  a'  for  some 
reciprocal  lattice  vector  a',*  a  condition  that  eq.  (33)  implies.  That  is,  the  condition  implied 
by  the  B&L  model  for  superradiance  in  a  crystal  is  sufficient  for  coupling  to  a  Borrmann 
mode,  the  propagation  direction  of  which  will  be  determined  by  the  ratio  of  the  wavelength 
X  to  the  atomic  spacing  d  in  accordance  with  eq.  (33). 


*  Then,  c.g.,  a  «  2a'  will  satisfy  (34). 
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m.  USE  OF  A  HEAT  SINK  TO  PROTECT  THE  LASING 
MATERIAL  AFTER  UPCONVERSION 


A.  INTRODUCTION 

It  has  been  recognized  for  a  long  time  that  heating  during  the  pumping  stage  would 
pose  a  serious  problem  for  a  y-ray  laser  (Ref.  27).  Avoiding  high  temperatures  is  crucial, 
since  the  Mdssbauer  Effect  will  be  needed  to  preserve  coherence  among  the  photons.  This 
restricts  the  temperatures  of  operation  to  well  below  the  Debye  temperature.  One  recent 
proposal  to  get  around  this  suggests  looking  for  an  isomer  with  a  nearby  state  at  slightly 
higher  eneigy  (Ref.  28).  The  isomer  can  be  slowly  pumped  without  heating  problems,  and 
the  lasing  then  triggered  by  rapid  "upconversion"  to  the  nearby  state.  This  would  lead  to  a 
single-pulse  laser. 

We  have  shown  previously  that  upconvCTSion  with  100  eV  to  10  keV  photons  will 
lead  to  catastrophic  heating  via  the  photoelectric  effect  (Ref.  29).  We  have  also  presented 
preliminary  calculations  which  indicate  that  the  heat  cannot  effectively  be  conducted  away 
(Ref.  30).  In  this  paper  we  show  that  the  heating  problem  gets  progressively  more  severe 
with  decreasing  photon  energy  all  the  way  down  to  the  ionization  threshhold. 
Furthermore,  in  this  chapter  and  the  next  we  confirm  that  conductive  cooling  of  the  lasing 
material  does  not  appear  to  provide  a  solution.  The  essence  of  the  problem  is  that  heat 
sinks  with  sufficient  heat  capacity  to  store  the  heat  from  the  photoelectric  heating  of  the 
lasing  material  are  so  large  they  generate  and  trap  enough  photoelectrons  themselves  to  heat 
up  to  temperatures  where  they  are  useless. 

Before  presenting  the  calculations  in  the  next  section,  the  physics  underlying  the 
results  will  be  sketched.  There  are  three  dimensionless  ratios  that  largely  control  heating 
effects.  The  first  is  the  ratio  of  the  energy  required  per  atom  to  melt  a  sample  or  heat  it  to 
the  Debye  temperature  compared  to  the  energy  of  the  nuclear  excitation  in  the 
upconversion.  Specific  heats  of  melting  are  orders  of  magnitude  smaller  than  the  lowest 
known  nuclear  transition  energy  (73  eV).  We  will  consider  energies  as  low  as  the 
inonizadcm  threshold,  but  these  energies  are  still  large  ccnnpared  to  the  energies  required  for 
melting  or  heating  to  the  Debye  temperature. 


The  second  and  third  ratios  characterize  the  relative  probability  of  upconversion  of 
the  nucleus  versus  photoionization  of  the  atom.  For  ideal  beams  tuned  to  the  nuclear 
resonance  energy,  the  relevant  ratio  would  be  the  photoelectric  cross  section  divided  by  the 
nuclear  resonant  cross  section.  For  ordinary  nuclear  transitions,  involving  keV  energies  or 
greater,  very  small  leakage  of  the  upconversion  energy  will  destructively  heat  the  sample 
consequently.  For  very  low  energy  transitions,  should  any  exist,  nuclear  resonance  cross 
sections  are  very  small,  and  the  photoelectric  cross  sections  comparable  or  even  larger. 
MOTe  energy  goes  into  photoelectric  heating  than  into  nuclear  upconversion. 

Finally,  in  the  real  world,  photon  beams  have  finite  widths.  These  beams  are  very 
broad  compared  to  the  nuclear  states.  The  ratio  of  the  photoelectric  excitation  to  nuclear 
excitation  will  always  be  higher  that  the  ratio  of  photoelectric  to  resonant  cro'^s  sections  by 
approximately  the  ratio  of  beam  width  to  resonance  width.  For  current  beams,  this  leads  to 
enormous  effective  increase  in  the  photoelectric  cross  section. 

Addition  of  a  heat  sink  to  conduct  away  the  energy  fails  because  the  heat  sink  itself 
is  heated  to  high  temperature  by  photoelectrons  produced  within  either  the  heat  sink  or  the 
lasing  material.  Basically,  viewed  as  a  function  of  the  radius  of  the  heat  sink,  by  the  time 
the  heat  sink  is  large  enough  to  handle  the  heat  and  electrons  from  the  lasing  material  and 
still  remain  at  low  temperature,  it  is  so  large  that  the  direct  photoelectric  heating  of  the  heat 
sink  is  fatal. 

B.  CALCULATIONS 

We  want  to  consider  the  temperature  rise  in  a  heat  sink  used  for  temperanue  control 
of  the  lasing  material  in  a  y-ray  laser.  The  essential  properties  of  the  heat  sink  are  good 
thermal  conductivity  and  low  photoelectric  cross  section,  with  the  latter  being  most 
important.  For  this  reason,  we  have  considered  lithium  as  our  heat  sink  material.  For  the 
lasing  material  we  continue  to  base  our  model  on  the  14.4  keV  Mossbauer  transition  in 
Fe57  (Ref.  29).  Note  the  calculations  use  the  solids  with  the  lowest  photoelectric  cross 
section  and  the  highest  Debye  temperature,  which  is  the  optimum  combination. 

In  this  section,  we  will  compute  the  temperature  of  a  Li  heat  sink  itself,  under 
various  assumptions  about  the  upconversion  transition  in  the  lasing  material.  There  are 
four  sources  of  heat.  First,  heat  produced  by  the  photoelectrons  in  the  iron  which  is 
conducted  to  the  heat  sink.  Second,  if  the  photon  energy  is  high  enough,  the  photo¬ 
electrons  produced  in  the  Fe  will  not  be  trapped  there.  They  will  continue  out  into  the  Li, 
where  they  will  continue  to  deposit  energy.  Assuming  no  change  in  dE/dx,  the  energy 


deposition  will  be  linear  in  the  path  length  in  the  lithium.  The  energy  deposition  per  atom, 
or  the  temperature  increase,  will  vary  as  the  reciprocal  of  the  diameter  of  the  heat  sink. 
These  first  two  contributions  favor  a  large  heat  sink.  The  direct  heating  due  to  photo¬ 
electrons  produced  within  the  lithium  itself  increases  quadratically  with  the  diameter  of  the 
heat  sink,  and  favors  a  small  heat  sink.  Finally,  as  the  residual  ion  from  which  the  electron 
was  removed  returns  to  its  gound  state,  energy  from  photons  and  Auger  electrons  is 
deposited.  These  contributions  are  assumed  to  have  radiation  lengths  characteristic  of 
100  eV  photons,  and  have  a  total  energy  deposition  limited  to  the  initial  photon  energy 
minus  the  photoelecmm  energy. 

The  baseline  calculation  assumed  a  0.0 Ip.  Fe  wire  embedded  in  a  Li  heat  sink. 
Borrmarm  reduction  of  the  photoelectric  effect  by  a  factor  of  100  was  assumed  for  both 
materials.  The  probability  of  upconversion  of  the  Fe^”^  was  held  constant  at  1  percent.  The 
basis  for  these  choices  is  described  in  Ref.  29. 

Assuming  the  number  of  lithium  atoms  is  large  compared  to  the  iron  atoms,  the  total 
energy  deposited  per  atom  of  Li,Q,  can  be  written 

,dEx  ^  Ppe  \ 

^  Pu  ^ 

^(Fe)ap^dx^Li  2 

+  (Ey-  Ee)  (l-exp(-aja^  )j  • 

PpE  is  the  probability  of  producing  a  photoelectron  per  iron  atom  in  the  lasing  sample.  For 
a  given  size  iron  san:q)le,  this  depends  on  the  beamwidth,  and  the  energy  of  the  transition. 
dE/dx)y  is  the  energy  lost  by  the  electron  per  unit  length  in  material  y  (Ref.  31). 
Technically,  the  two  values  for  Li  should  be  different,  since  iron  and  lithium  produce 
different  energy  photoelectrons.  The  smaller  value  was  used  in  all  cases  in  keeping  with 
our  neglect  of  the  increase  in  dE/dx  as  the  electron  slows  down.  The  radius  or  width  of  the 
material  y  is  ay,  and  the  number  density  py  the  reciprocal  mean  free  paths  for  the  photo¬ 
electrons  labeled  ii(y),  and  the  radiation  length  characterizing  the  low-energy  photons  and 
Auger  electrons  is  urad  (Ref.  31). 
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In  the  above  expression,  the  four  terms  correspond  to  the  four  physical  sources 
discussed  above.  The  optimum  value  of  a  can  be  found  by  differentiation.  For  incident 
photons  with  about  10  keV  energy,  the  second  and  third  terms  are  dominant.  This  leads  to 
the  approximate  result 


It(Li)  pLi 


=  0.32p  . 


A  numerical  solution  to  the  minimum  heating  leads  to  values  of  aLi/ape  around  30, 
so  we  used  au  =  0.30  pm  in  these  calculations.  The  overall  energy  scale  is  set  by  the 
probability  of  the  photoelectric  effect  occuring  per  iron  atom.  This  depends  on  the  size  of 
the  iron  sample,  the  upconversion  percentage,  the  relative  sizes  of  upconversion  and 
photoelectric  cross  sections,  the  Borrmann  reduction  of  the  photoelectric  effect,  and  the 
photon  beamwidth.  We  assume  a  0.01  pm  Fe  sample,  and  Borrmaim  reductions  of  100  in 
both  Fe  and  li.  These  choices  are  discussed  in  our  earlier  work  (Ref.  29).  The  physical 
cross  sections  are  known  (Ref.  32).  The  essential  remaining  parameter  is  the  beam  spread, 
or  more  properly,  the  ratio  of  the  beam  spread  to  the  width  of  the  nuclear  state.  We  have 
used  10^  as  our  baseline  value.  By  comparison,  the  width  of  the  14  keV  state  in  Fe57  is 
seven  orders  of  magnimde  smaller  than  the  smallest  beam  spread  ever  achieved  in  this 
energy  region. 

Using  10^  as  our  value  for  the  ratio  of  beam  ^read  to  intrinsic  state  width,  we  find 
approximately  0.25  eV/atom  energy  depositon  in  Li,  corresponding  to  a  temperature  of 
around  8(X)  K.  This  teii^>erature  is  well  above  the  Debye  temperature  of  any  material,  and 
clearly  does  not  represent  a  solution  to  the  heating  problem.  Improved  beam  resolution 
does  allow  for  a  solution.  For  a  beam  spread  of  10^  times  the  intrinsic  state  width,  the 
heating  would  be  minimal,  and  the  heat  sink  would  work.  (See  more  detailed  calcualtion  of 
this  case  in  the  next  section.)  However,  the  improvements  in  beam  resolution  must  be 
accompanied  by  dramatic  increases  in  intensity  before  this  becomes  feasible.  Current 
beams  have  only  a  few  tens  of  photons  per  typical  nuclear  state  width. 

As  the  energy  decreases,  the  photoelectric  cross  sections  increase,  except  at  shell 
crossings.  The  rate  of  energy  loss  increases,  and  the  width  of  the  nuclear  state  decreases. 
These  features  all  tend  to  make  the  heating  get  worse  as  the  energy  decreases.  For  the  case 
under  consideration,  the  peak  resonance  and  photoelectric  cross  sections  are  plotted  as  a 
function  of  the  resonance  energy  in  Fig.  6.  The  peak  resonance  cross  section  is  computed 
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assuming  a  nuclear  structure  matrix  element  equal  to  that  of  the  14  keV  transition  in  Fe^^, 
and  incoiporating  the  energy  cubed  phase  space  factor  in  the  excitation  width  and  the 
wavelength  squared  factor  in  the  cross  section. 
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Figure  6.  Resonance  and  Photoelectric  Cross  Sections 


With  the  exception  of  the  shell  crossings,  the  heating  problem  gets  progressively 
worse  as  the  energy  decreases.  Even  at  the  K-shell  crossing,  the  situation  is  no  better  than 
at  10  keV.  The  next  physical  threshold  at  which  this  situation  could  change  is  the 
ionization  threshold,  at  around  10  eV.  Below  that  energy,  it  might  be  possible  to  upconvert 
a  nuclear  state  with  a  laser  if  such  a  low-energy  nuclear  transititm  existed. 

C.  CONCLUSIONS 

Upconversion  of  an  isomer  to  a  nearby  excited  state  for  the  purposes  of  initiating 
y-ray  lasing  has  been  investigated  from  the  point  of  view  of  photoelectric  heating  of  the 
lasing  material.  This  has  long  been  known  to  be  a  serious  problem,  because  it  is  necessary 
to  maintain  temperatures  low  compared  to  the  Debye  temperature  in  order  to  exploit  the 
Mbssbauer  Effect  and  achieve  coherence. 
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In  the  10  keV  regime,  nuclear  resonance  cross  sections  are  large  compared  to 
photoelectric  cross  sections.  Nevertheless,  available  beams  are  so  wide  compared  to 
nuclear  intrinsic  state  widths  that  the  photoelectric  heating  would  destroy  a  sample  before 
1  percent  upconv^sion  was  possible.  This  is  true  even  with  substantial  suppression  of  the 
photoelectric  cross  sections  via  the  Borrmann  Effect.  To  achieve  upconversion  without 
excess  heating  would  require  a  beam  with  10^^  photons/burst  with  a  0.(X)5  eV  spread.  As 
the  energy  decreases,  the  requirements  on  the  spectral  resolution  become  more  severe,  until 
the  ionization  threshold  is  crossed  and  the  photoelectric  effect  vanishes. 

Upconversion  of  nuclear  states  in  the  10  eV  to  10  keV  energy  regime  will  not  be 
possible  without  improvements  of  several  orders  of  magnitude  in  both  beam  intensity  and 
beam  resolution.  This  conclusion  caimot  be  evaded  by  resort  to  a  heat  sink,  because  any 
heat  sink  large  enough  to  remove  sufficient  heat  from  the  lasing  material  would  be  so  large 
that  it  traps  too  much  heat  itself. 


IV.  BLACKBODY  NUCLEAR  SPECTROSCOPY 


A.  INTRODUCTION 

One  can  inomgine  two  approaches  to  building  a  y-tay  laser.  One  could  induce  decay 
by  manipulating  the  lifetime  of  a  long-lived  state,  or  one  could  pump  nuclei  up  to  a  shon- 
lived  state.  One  difficulty  with  the  second  process  is  that  power  requirements,  and 
attendant  heating  and  sample  destruction  problems,  make  pumping  of  nuclear  levels 
infeasible.  This  can  be  partly  evaded  by  using  an  isomer,  which  can  be  pumped  slowly, 
widiout  excessive  heating  or  power  requirements,  in  order  to  store  most  of  the  energy  to  be 
released  by  the  laser.  Subsequently,  lasing  can  be  initiated  by  another  pumping  step, 
sometimes  called  upconversion,  which  is  easier  to  accomplish. 

When  upconversion  firom  an  isomer  is  examined  mtne  closely,  it  becomes  clear 
that,  at  least  for  transitions  in  the  10  eV  to  10  keV  region,  heating  of  the  sample  due  to  the 
photoelectric  effect  is  still  aproblem.  We  believe  this  is  a  solid  result,  demonstrated  by  botii 
"back-of-the-envelope"  and  detailed  analyses  repotted  previously.  Nevertheless,  the  result 
is  controversial;  we  believe  this  because  other  workers  have  asserted  that  in  the  SO  keV 
region  heating  will  not  be  a  problem.  We  suspect  that  the  heating  problem  at  higher  x-ray 
enerpes  will  depend  strongly  on  Z,  and  may  hinge  on  the  feasibility  of  conductive  cooling 
witii  low  Z  materials.  These  are  worthwhile  questions  for  furtiier  analysis. 

Here,  however,'  we  look  at  the  other  end  of  the  energy  spectrum.  Below  about 
10  eV  the  photoelectric  cross  section  vanishes  because  the  available  energy  is  below  the 
ionization  potential.  As  a  result,  the  heating  will  drop  precipitously.  This  has  motivated  us 
to  consider  how  one  would  find  nuclear  levels  with  short  lifetimes  located  near  isomers. 
By  "near"  in  this  context  we  mean  below  the  ionization  threshold,  so  we  are  talking  about 
"optical  nuclear  transitions." 

The  existence  of  nuclear  energy  levels  with  short  lifetimes  lying  near  longer-lived 
isomeric  levels  has  not  yet  been  shown  experimentally,  due  to  the  present  limitations  of 
nuclear  spectroscopy,  but  theoretical  models  anticipate  their  existence.  These  short-lived 
states  are  essoitial  to  proposed  two-stage  pumping  schemes  for  the  development  of  y-ray 
lasers,  where  the  nucleus  is  "pumped"  from  the  storage  isomer  to  the  nearby  level  from 
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which  the  lasing  transition  can  take  place.  For  this  reason  and  for  pure  research  purposes, 
various  studies  have  been  done  cm  possible  methods  of  locating  such  energy  levels. 

Previous  studies  on  the  feasibility  of  voifying  die  existence  of  these  transfer  levels 
have  included  work  done  by  R.C.  Haight  and  G.C.  Baldwin  (Ref.  33)  of  the  Los  Alamos 
National  Laboratory  (LANL),  and  EC.  Zimm»mann  (Ref.  34)  of  the  Institute  for  Defense 
Analyses  (IDA).  The  Baldwin  and  Haight  study  looked  at  the  prospect  of  using  broadband 
radiation  to  excite  the  transfer  level.  Short  pulses  of  blackbody  radiation  produced  by  a 
laser-target  interaction  would  be  used  to  illuminate  the  isomer,  and  a  small  band  of  this 
spectrum  would  convert  the  isomer  to  the  transition  level.  It  was  concluded  by  Haight  and 
Baldwin  that  this  mediod  would  prove  to  be  too  restrictive  to  be  widely  useful.  In  order  for 
a  sufficient  amount  of  the  isomer  to  be  converted  to  the  transfer  level  in  10~^  s  (the  length 
of  typical  laser  pulses),  the  multipolarity  of  tiie  transition  could  only  be  El  or  Ml,  the 
transition  strengths  would  need  to  be  close  to  the  single  particle  value,  and  the  transition 
energy  would  be  limited  to  the  range  greater  than  1(X)  eV  for  the  transition  strengths  to  be 
laiige  enough  and  less  than  10  keT  for  the  phoum  density  to  be  high  enough. 

The  method  proposed  by  Zimmermann  requires  heating  long-lived  Mtissbauer 
isotopes  to  high  temperatures.  This  should  thermally  excite  a  sufficient  amount  of  the 
isomer  to  the  transfer  level  to  noticeably  increase  the  rate  of  decay  m  ground  state.  An 
estimate  for  the  value  of  AE  (the  energy  separation  between  the  two  levels)  could  be  found 
by  slowly  increasing  tite  temperature  from  3(X)K  and  noting  increase  in  the  decay  rate  witii 
temperature.  This  mediod  is  based  oa  the  assunqition  that  rtqnd  diermal  equilibration  of  the 
two  excited  states  holds,  and  if  this  assun^tion  is  inconect,  the  increase  in  the  decay  rate  at 
elevated  temperatures  might  not  be  noticeable. 

This  study  examines  a  potential  method  fiv  locating  transfer  levels  which  uses  a 
cavity  (either  a  blackbody  cavity  ot  a  coherent  source  cavity)  filled  with  a  dilute  atomic 
vapor,  containing  the  requisite  isomeric  energy  levels.  A  vaporous  substance  is  used  so  as 
to  avoid  unwanted  solid-state  material  effects.  The  isomers  will  be  illuminated  with  either 
the  continuous  photon  plectrum  associated  with  the  radiation  fiom  a  blackbody,  or  with  the 
radiation  produced  by  a  ctdierent,  tunable  source  (such  as  a  microwave  power  source  or  a 
laser).  Either  of  these  methods  of  photoexcitation  should  produce  an  intense  photon  flux, 
and  botii  methods  would  be  relatively  chetp  and  easy  to  use. 

For  ease  of  calculation,  this  study  has  chosen  a  doubled  count  rate  as  a  benchmark 
value,  i.e.,  a  short-lived  nearby  energy  level  is  assumed  to  have  been  successfully  located 
when  die  induced  decay  rate  is  equal  to  the  spontaneous  decay  rate. 
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This  method,  unlike  that  proposed  by  Haight  and  Baldwin,  allows  for  a  longer 
count  time  and  therefore  a  oxne  gradual  heating  of  the  isomer.  This  study  assumes  the 
isomers  will  be  heated  to  between  1000  K  and  5000  K.  Unlike  the  method  proposed  by 
Zimmermann,  the  thermal  equilibration  is  calculated  rather  than  assumed.  Another  possible 
advantage  to  using  this  method  is  that  the  energy  of  the  transfer  state  need  not  be  known  to 
be  within  any  specific  energy  range,  it  can  be  searched  for  over  a  tnoad  range  of  values. 
Transfer  levels  for  more  than  one  isomeric  level  can  be  searched  for  simultaneously,  as 
long  as  the  energies  of  the  isomers  associated  with  the  various  target  nuclei  are  distinct 


For  a  blackbody  cavity,  the  count  rate  is  given  by 


dk  , 


where  cn^  is  the  photon  flux  in  photons/cm^-s-Ak,  which  for  a  blackbody  is 


cn  _ x  =  ^  =  -^ 

,2  e«-l  •  k,T  kjT 


and  0(k)  is  the  Breit-Vi^gner  cross  section  for  low'energy  photon  scattering 


o(k) 


no 

2n 


r.r. 


where  Fe,  Fi,  and  Ft  are  die  elastic,  inelastic,  and  total  widths;  E  and  Eq  are  the  energies  of 
the  particle  and  the  resonance;  and  X  and  Xo  are  the  wavelengths  of  the  particle  and  the 
resonance.  Therefore,  the  count  rate  is 


dN^fk^  c  ^^eF;  ^ 

*  J  2ir^  c*-l  ^_E^)2  +  (F/2)^ 


We  assume  that  the  inelastic  width  is  approximately  equal  to  the  total  width  (F  =  Ft),  and 
that  the  elastic  width  is 

FeSZlxio'^E^  , 


where  By  is  the  decay  energy  in  eV.  This  assumes  an  Ml  transition  with  a  strength  of 
1  Weisskopf  Unit,  characteristic  of  the  expected  width  of  a  strong  single-particle  magnetic 
dipole  transiticm. 

The  count  rate  integral  will  have  non-negligible  contributions  at  k  =  0  (E  =  0)  and  at 
resonance  (E  »  Eq).  This  motivates  us  to  divide  the  total  contribution  into  two  pieces,  i.e.. 
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.  The  motivation  for  this  is  that  the  low-frequency  and  on- 


dN  _  dNres  dNiao 
dt  “  dt  dt 

resonance  contributions  are  physically  quite  different,  and  the  low-frequency  modes  iiught 
be  eliminated  or  damped  in  a  plasma.  If  possible,  we  would  want  to  neglect  them.  This 
proves  to  be  sensible  in  the  region  4kT  <  Eq  <  20  kT;  for  E©  above  4  kT  there  is  a 
minimum  between  the  zero-energy  and  resonance  "peaks."  For  E©  above  20  kT,  the 
exponential  decrease  of  tik  begins  to  overwhelm  the  resonance  contribution. 


Assuming  that  (E  -  Eq)  »  Ft,  (i.e.,  far  off  resonance),  the  zero-energy  peak 
ccxitiibution  to  the  count  rate  is  iq>pioximately 

,2  .  .  ^-20, 


=^(kB T)^ (2.1  X  10 


This  result  is  obtained  by  replacing  the  full  expression  for  0(k)  with  the  zero-energy 
coefficient  of  1/k,  and  evaluating  the  resulting  integral  analytically.  This  eliminates  the 
resonance  peak  at  E  =  Eo  from  the  calculation.  Since  the  enhanced  decay  rate  is  inversely 
proportional  to  the  mean  lifetime  of  the  excited  state,  in  order  to  double  the  spontaneous 
decay  rate  we  would  need 

(2.2x10-^)  , 

where  1/x*  is  the  reciprocal  of  the  excited  state  lifetime,  which,  according  to  the  energy¬ 
time  uncertainty  principle,  satisfies  the  relation  l/v  ^  F/h.  Therefore,  we  have 

_ af(6.ixio")  (looon-)^  • 


Thus,  we  would  require  an  isomer  with  a  mean  lifetime  some  20  orders  of  magnitude 
greater  than  the  excited  state  lifetime.  In  general,  the  greatest  increase  in  decay  rate  that 
could  be  expected  would  occur  if  the  Ml  transidcm  connected  nearby  states  where  the  lower 
state  had  a  primary  decay  mode  of  multipolarity  M  (L  +  1)  and  the  upper  state  decay  by 
multipolarity  E  (L)  (see  sketch  below). 
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Again,  using  the  Weisskopf  estimates  we  find  for  the  ratio  of  radiative  widths 
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if  both  transitions  are  of  energy  Ey  (in  eV).  In  order  to  achieve  the  needed  enhancement  we 
require 
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which  means  that  the  transition  energy  of  the  isomer  to  the  ground  state,  Ey,  must  be 

Ey^ 


^  ~  £  0.1  eV  for  T  =  2500  K  . 
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Grouped  states  would  thus  be  required  in  the  sense  that  the  isomer  and  the  higher-lying, 
rapidly-decaying  state  would  both  need  to  be  very  close  to  a  third  state  to  which  they  both 
could  decay.  This  appears  sufficiently  implausibie  to  us  that  we  will  not  pursue  it  further. 

The  contribution  of  the  decay  rate  integral  from  near-resonance  (within  m  widths)  is 

^■MnT 

dN„.  f 

Etf-tnT 


Provided  only  that  ksT  is  large  compared  to  the  total  width,  n^  can  be  factored  from  the 
expression.  The  count  rate  can  then  be  approximated  by  the  full  integral 

^  =  cn^(10  Jo(k)dk.4(2.1xlO-“)i-|;;S_  . 

—  e  ^  “1 

Figure  7  shows  the  lifetime  (reciprocal  of  the  count  rate)  as  a  function  of  the  energy  for 
values  of  the  temperature  ranging  from  1000  K  to  5000  K.  The  increase  in  lifetime  at 
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Figure  7.  Lifetime  ve.  Energy  for  Varioue  Temperatures 
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lower  energies  is  due  to  the  decrease  in  the  Breit-Wigner  cross  section  which  goes  to  zero 
as  the  energy  of  the  transition  goes  to  zero.  The  increase  at  higher  energies  is  due  to  the 
decreasing  photon  flux,  which  is  a  characteristic  of  the  blackbody  curve  at  high  energies. 
For  a  blackbody  source  at  5000  K,  an  excited  state  lifetime  of  1  year  corresponds  to  an 
energy  of  approximately  6  eV.  Figure  8  shows  the  lifetime  as  a  function  of  the  energyfor  a 
blackbody  at  5000  K  and  a  blackbody  at  50,000  K.  This  curve  supports  the  Baldwin- 
Haight  conclusion  that  the  use  of  short  pulses  of  blackbody  radiation  would  not  be  useful 
in  searching  for  the  desired  short-lived  energy  levels.  A  lifetime  of  10“^  s  would  be  well 
off  the  bottom  of  the  curve  shown,  requiring  an  unreasonably  high  blackbody  temperature. 
From  this  we  conclude  that  photoexcitation  using  a  blackbody  source  would  be  a  useful 
method  for  searching  for  low-excitation  transitions,  based  on  an  isomer,  for  energies  from 
0.1  to  10  eV. 

The  results  of  tiiese  blackbody  calculatirais  can  be  summarized  as  follows: 

1 .  The  contribution  to  the  excitatirai  is  dominated  by  the  photon's  near-resonance. 

2.  A  nuclear  state  approximately  0.1  eV  to  10  eV  above  an  isomer  could  be 
detected  in  a  blackbody  cavity  if  it  had  a  fairly  strong  Ml  transition  connecting 
it  to  the  isomer.  This  depends  in  detail  on  the  strength  of  the  transition  and  the 
lifetime  of  die  isomer. 

3.  These  results  are  consistent  with  tlK>se  of  Baldwin  and  Haight  in  the  sense  that 
even  at  50,000  K,  there  is  scant  probability  of  a  count  after  1  ns. 

For  photoexcitation  using  a  coherent  power  source  (e.g.,  a  microwave  power 
source),  we  assume  that  the  photon  flux  follows  a  Gaussian  distribution  such  that 

7  ^(k-k„yw/ 

n^  =  CjJ  e  dk  , 

where  Ci  is  a  normalization  constant  and  is  the  width  of  the  peak.  Integrating  over  all 
possible  wave  numbers,  we  Bnd  that  the  total  photon  flux  is 

n^  =  CjVnW^  . 

The  energy  flux  rate  is  the  product  of  the  photon  flux  and  the  energy  per  photon: 

energy  flux  rate  =  Cj  Vn  Wj^Ticko  . 
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Figure  8.  Lifetime  vs.  Energy  for  Temperstures  of  5000  a;  50,000  K 
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Equating  this  with  the  power  per  cross-sectional  area  of  a  microwave  cavity,  Ci  is  found 
to  be 


(  'I 

1.60  X  10"^’  J 

V  y 

yjc  Wj^h ckj, 

^  > 

where  P  is  the  power  of  the  source  in  watts,  and  A  is  the  cross-sectional  area  of  the  source 
cavity.  For  a  1-kW  power  source  with  a  cross-sectional  area  of  100  cm^,  Ci  is 

27  2 

approximately  1.78  x  10  /k^ .  One  of  the  factors  of  ko  comes  from  the  assumption  that 
the  width  of  the  peak  is  a  function  of  ko  (we  used  =  ko/10(X))  and  the  other  factor  of  ko 
comes  from  the  energy  per  photon,  E  =  hcko.  Therefore,  the  total  induced  decay  rate 
corresponding  to  the  use  of  a  coherent  power  source  is 

^  =  (5.24  X  10^)  c  h^  /  Eo  . 
dt 

The  possibility  of  using  other  types  of  coherent  sources,  such  as  lasers,  should  also 
be  considered.  Figure  9  shows  the  lifetime  as  a  function  of  energy  for  various  types  of 
sources.  Included  are  a  1-kW  coherent  power  source,  characteristic  of  a  commercially 
available  microwave,  and  various  types  of  continuous  wave  lasers,  along  with  the  ranges 
of  energies  over  which  they  are  tunable. 

We  find  that  below  about  1  eV  a  tunable,  relatively  high  energy  resolution  (i.e., 
coherent)  source  would  be  more  effective  than  using  blackbody  radiation.  At  very  low 
energies,  microwave  sources  would  suffice.  Near  1  eV,  IR  lasers  would  be  needed.  For 
this  search  not  to  be  hopelessly  tedious,  it  would  be  necessary  that  lasers  be  available 
which  were  tunable  over  a  few  percent  of  their  centroid  energy,  at  least  At  present  tunable 
lasers  are  "tunable"  over  a  few  natural  line  widths.  This  situation  is  expected  to  improve, 
perhaps  in  the  next  five  years. 

B.  SUMMARY 

Upconversion  schemes  for  a  y-ray  laser  would  be  greatly  aided  by  near-degeneracy 
of  nuclear  levels  connected  by  an  Ml  transition  of  appreciable  strength.  We  refer  to  such  a 
transition  as  an  optical  nuclear  transition.  We  have  investigated  mechanisms  for  locating 
such  states.  At  present,  the  best  technique  would  be  to  heat  a  sample  of  the  isomer  in  a 
blackbody  cavity  and  look  for  enhanced  decays  as  a  result  of  photoexcitation  of  the  nearby 
state.  This  would  be  feasible  in  the  energy  region  from  0.1  eV  to  10.0  eV. 
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For  energies  below  about  0.1  eV,  the  cross  section  is  too  small  for  this  approach  to 
be  feasible.  The  transition  could  be  driven  with  commercially  available  coherent  sources,  if 
the  energy  was  right.  Coherent  sources  tunable  over  a  wide  range  are  expected  in  the  fairly 
near  future,  at  which  time  the  search  for  optical  nuclear  transitions  could  proceed  down  to 
very  low  energy. 

The  efforts  that  would  be  required  using  coherent  sources  at  present  seem 
unwarranted.  Driving  the  transitions  with  blackbody  radiation,  however,  should  be  a  fairly 
cheap,  easy  experiment  which  might  produce  some  interesting  nuclear  spectroscopy,  in 
addition  to  advancing  the  feasibility  of  the  y-ray  laser. 
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V.  ENERGY  DISSIPATION  DURING  PUMPING  OF 

y-ray  laser 


A.  INTRODUCTION 

Upconversion  of  an  isomeric  nuclear  level  to  a  short-lived  nearby  level  suitable  for 
stimulated  emission  is  currently  the  most  popular  concept  for  producing  inversion  in  a 
nuclear  system  (Ref.  3).  The  pumping  mechanism  can  be  either  accomplished  through  a 
direct  interaction  of  electromagnetic  radiation  with  the  nucleus  (Ref.  35)  or  through  an 
intermediate  step  involving  first  laser  exdtadon  of  electrons  followed  by  energy  transfer 
between  atomic  electrons  and  the  nucleus  (Ref.  36).  For  both  mechanisms,  large  electronic 
absorption  of  the  laser  beam  could  produce  material  modification  resulting  at  least  in  a 
reduction  of  the  Mdssbauer  and  Borrmann  effects  and  possibly  in  destruction  of  the 
crystalline  structure  of  die  sample.  Sparrow  (Ref.  37)  has  estimated  the  heating  effects  as  a 
function  of  punqiing  beam  intensity,  width,  and  energy  using  thermodynamic  arguments. 
He  also  calculated  the  effect  of  a  heat  sink  in  reducing  the  temperature  rise  (Ref.  38). 

The  main  purpose  of  this  chapter  is  to  extend  the  previous  calculations  of 
D.  Sparrow  (Refs.  30,  37,  38)  to  more  realistic  geometries,  including  different  material 
combinatitms  and,  most  importantly,  to  calculate  dynamic  or  transient  effects  and  study  the 
approach  to  steady  state.  The  second  purpose  is  to  use  a  general  approach,  applicable  by  a 
change  of  one  input  parameter,  to  rectangular,  cylindrical  and  spherical  geometries,  to  the 
analysis  of  heat  ccmduction  problems  specialized  to  the  cylindrical  case. 

This  approach  is  valid  for  general,  time-varying,  extended  heat  sources.  This 
capability  will  be  useful  to  investigate  the  heating  problem  and  the  effectiveness  of 
proposed  solutions. 

To  carry  out  this  program  we  address  the  general  problem  of  heat  generation  in 
specific  regims  of  a  system  and  its  transfer  to  other  areas  of  the  system  through  the  process 
of  diffusion.  The  computer  program  used  to  perform  the  calculations  is  based  on  a  finite 
element  technique  described  in  Appendix  D.  The  code  is  written  assuming  axial  symmetry 
but  otherwise  permits  rather  general  conditions.  The  heat  generation  and  thermal  properties 
of  die  media  can  vary  widi  time  and  space  in  an  arbitrary  way. 
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Consider  an  arbitrary  geometry  with  the  medium  of  interest  composed  of  two 
regions.  One  region  (I,  the  active  volume)  is  composed  of  a  good  Mdssbauer  isotope  that 
absorbs  the  radiation  preferentially  in  a  nuclear  transition  with  some  radiation  absorbed  by 
the  atoms  leading  to  production  of  heat  in  the  sample.  Another  region  (11)  acts  as  a  heat 
sink  for  the  active  region  because  of  its  high  thermal  conductivity,  low  absorption  in  the 
wavelength  of  the  pumping  radiation,  and  high  specific  heat.  How  effective  can  such  a 
heat  sink  be  in  preventing  the  temperature  rise  in  the  active  volume  from  melting  the 
crystal?  This  question  has  been  addressed  systematically  by  Dr.  David  A.  Sparrow  who 
calculated  heat  generated  in  the  active  region  for  a  hypothetical  isotope  having  the  best 
Mdssbauer  properties  known  (Fe^^),  together  with  an  assumed  isomeric  storage  level 
within  eV  to  keV  of  the  Mdssbauer  upper  state.  Such  an  isomer  is  not  known  to  exist  and 
represents  an  extreme  optimistic  case.  According  to  other  calculations  (Sparrow,  Ref.  38), 
which  assumed  an  infinite  conductivity  of  the  sink  medium  (11),  the  heat  generated  would 
still  destroy  the  crystal.  It  was  decided  that  more  realistic  cases  should  be  treated  to 
examine  the  phenomena.  In  the  present  study,  we  performed  calculations  using  real 
thermal  properties  of  the  media  and  realistic  space  variations  of  the  heat  generation  in  the 
media  (I  &  11)  to  determine  under  what  ccmditicms  a  passive  heat  sink  can  protect  the  active 
region  from  being  destroyed  during  the  pumping  process. 

B.  CALCULATIONS/RESULTS 

1 .  Specific  Geometry 

The  calculations  were  performed  for  a  cylinder  of  one  material  (Fe^^)  encased  in  a 
sleeve  of  another  material  (Li).  The  geometry  is  shown  in  Fig.  10.  In  (a)  the  cylinder  of 
diameter  D]  is  the  active  region  and  the  sleeve  of  outer  diameter  D2  is  the  heat  sink.  In  (b) 
a  cut  containing  the  axis  of  the  cylinder  (shaded)  is  shown  with  the  four  test  points  at  which 
the  temperature  as  a  function  of  time  is  calculated  are  indicated.  In  (c)  two  assumed  heat- 
generation  profiles  which  are  functions  of  the  radius  x  are  shown. 

2.  Types  of  Calculations 

A  set  of  calculaticms  were  performed  with  an  initial  step  function  temperature  profile 
and  no  additional  heat  input  Generally,  the  temperature  was  assumed  to  be  100  K  initially 
in  medium  I  and  0  in  medium  11.  Both  temperature  profiles  as  a  function  of  radial 
direction,  and  temperature  changes  at  the  four  points  shown  in  Fig.  KKb)  are  graphed  in 
Figs.  11  and  12. 
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Figure  10.  System  Definition,  (a)  Geometry  used  in  the  caicuiations, 
(b)  Poaition  in  the  system  of  special  pointa  for  which  temperature 
as  a  function  of  time  is  plotted  in  subsequent  figures, 

(c)  Spatial  distribution  of  heat  sources. 
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A  set  of  calculations  were  perfonned  with  initial  teiiq}erature  0  °C  through  the 
medium  and  a  heat  input  (due  to  radiation  ionization)  assumed  in  each  medium. 
Two  different  radial  functions  of  the  heat  input  were  assumed  as  shown  graphically  in 
Fig.  10(c): 

(a)  1.03  X  10^  Cal/cm3/s  Region  I 
1.88  X  10“^  Clal/cm^/s  Region  II 

(b)  1.03  X  109  Cal/cmVs  Region  I 

7 

1,29^  10_  Cal/cm^/s  Region  n. 

3.  Calculated  Specific  Results 

Our  first  study  dealt  with  the  dynamics  of  the  approach  to  equilibration.  For  the 
system  shown  in  Fig.  10  with  Di  -  0.2  p.  and  D2  =  2.0  p  and  either  all  iron  or  iron  and 
lithium  combination.  Figs.  ll(a,b,c,d)  show  that  the  system  (inside  case),  r  <  Di/2 
tezzq)etature  decreases  to  about  10  percent  of  its  initial  ten:q)erature  in  just  over  1  ns.  Thus, 
heat  removal  by  passive  diffusion  effectively  occurs  on  the  time  scale  of  1  ns  for  this 
combination  of  geometry  and  thermodynamic  properties.  Any  attempts  at  active  heat 
removal,  such  as  flowing  coolants  past  the  whisker  combination  or  thermoelectric  cooling 
would  have  to  do  much  better  than  this  in  ordo^  to  be  useful 

When  we  introduce  a  distributed  heat  source  into  the  problem  there  is  an  initial  fast 
rise  in  temperature  followed  by  a  steady  state  situation  during  the  course  of  the  heating 
(pumping).  The  nonlinear  region  lasts  for  about  1  to  5  ns,  depending  on  the  exact 
conditions.  The  approach  to  equilibrium  is  exhibited  clearly  in  Fig.  12(k),  where  the 
temperature  differaice  in  time  at  a  particular  pcwt  is  plotted  as  a  function  of  time. 

Our  results  are  summarized  in  Fig.  13,  which  compares  the  steady  state  part  of  the 
heating  curves  for  (a)  an  iron  whisker  Di  =  0.1  p  and  no  sink  (D2  =  0)  and  two  heating 
curves  calculated  for  an  iron  whisker  with  a  lithium  sink  with  a  constant  heat  source  in  each 
region  (b)  and  a  radial  varying  heat  source  in  region  n  (c).  The  final  temperatures  after  a 
1  ps  pulse  are  513  ®C  for  (a),  76.0  for  (b),  and  64.2  ®C  for  (c). 

All  the  calculations  were  done  assuming  a  source  beamwidth  to  linewidth  ratio  of 
AE/T  s  10^,  which  is  a  very  optimistic  assumption.  Perhaps  a  nxne  realistic  assumption 
would  be  AE/T  *  10^,  which  would  increase  the  temperature  rise  by  two  orders  of 
magnitude. 
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C.  CONCLUSIONS 


A  heat  conduction  program  written  for  cylindrical  geometry  (axial  symmetry)  and 
capable  of  dealing  with  arbitrary  material  combinations  and  heat  source  locations  was  used 
to  examine  the  heating  problem  resulting  from  photon  pumping  of  a  y-ray  laser. 

Our  calculations  showed  that  for  a  system  on  the  order  of  microns  in  diameter, 
passive  heat  transfer  from  the  active  core  to  the  reservoir  occurred  on  the  order  of 
nanoseconds.  Any  active  technique  for  cooling  the  active  core  region  of  the  y-ray  laser 
would  have  to  be  faster  to  be  useful 

It  was  also  determined  that  after  a  few  nanoseconds  the  system  approached  a  steady 
heat  flow  situation  during  constant  heat  input  Where  comparisons  could  be  made  our 
results  connate  well  with  the  steady-state  idealized  geometry  calculations  of  Sparrow  (Ref. 
37  and  Chapter  IV  of  this  report). 


Figure  11.  (a)  Temperatures  at  points  1,  2,  3,  4  (positions  of  points  shown  In 
Fig.  10)  after  an  Initial  condition  of  100  °C  In  region  I  composed  of  Iron 
(Di  X  0.2  p)  and  0  °C  in  region  II  composed  of  IKhlum  (D2  s  2.0  p).  There 
are  no  heat  sources  In  I  or  II.  The  temperature  profiles  at  several  times  for 
this  case  are  shown  In  (b)  (when  the  thickness  of  region  II  Is  doubled  the 
results  do  not  change  substantially). 

(c)  The  same  situation  as  In  (a),  except  region  II  Is  composed  of  Iron. 

(d)  Temperature  profiles  for  (c). 
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Figure  12.  (a)  Temperature  changes  at  points  1,  2,  3,  4  assuming  Initial  T  s  0 

In  both  regions  and  a  heat  Input  of  1.03  x  10^  Cal/cm^/s  In  region  I  (Iron), 

and  1.88  X  10^  Cal/cm^/s  In  region  il  (lithium)  are  presented. 

(b)  Temperature  proflIes  for  selected  times  of  case  (a)  are  shown. 

(c)  Same  as  In  (a)  except  the  heat  input  in  region  ll  Is  a  function  of  radius, 

1.29  x10*/Cal/cm^» 

X  (in  cm) 

(d)  Temperature  profiles  for  selected  times  of  case  (b). 

(e)  Same  as  (c)  except  the  conductivity  in  region  II  Is  increased  by  a  factor 
of  five  (hypothetical  material  Indicated  as  LI*)  without  changing  the 
other  properties. 

(f)  Temperature  profile  for  selected  times  of  case  (e). 

(g)  Same  as  (c)  except  the  specific  heat  In  region  11  Increased  by  a  factor 
01  five  (hypothetical  material  Indicated  as  LI**)  without  changes  in  the 
other  properties. 

(h)  Temperature  profiles  for  selected  times  of  case  (g). 

(I)  Same  as  (c)  except  heat  Input  In  region  II  is  Increased  by  an  order  of 
magnitude  because  region  II  starts  at  x  s  o.Ol  Instead  of  O.i  and 

1.29  x10*/Cel/cni*/s 
X  (In  cm) 

(J)  Temperature  profiles  for  case  (h)  at  selected  times. 

(k)  Change  In  temperature  at  points  (1)  as  a  function  of  time 

[At  s  T(t  +  At)  -  T(t)]  for  cases  12(a)  and  12(b),  showing  the  approach  to 
steady  state. 
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Temperature  in  at  center  of  cyinder 


Figure  13.  (Summary  figure)  Temperature  change  in  central  iron  cylinder  (pt.  1) 
under  different  conditions,  assuming  the  initiai  temperature  is  T  s  0  for 
both  regions  and  Df  =  0.2  p,  D2  s  2.0  p.  The  solid  curve  shows  the 
temperature  when  the  heat  inputs  are: 

1.03  X  10^  Cai/cm^/s  in  region  I,  and 
1.88  X  10^  Cal/cm^/s  in  region  II 

The  dashed  curve  gives  the  temperatures  when  the  heat  input  is  a  function 
of  radius  in  region  II,  i.e.,  1.29  x  lO^/Cal/cm^/s.  The  dotted  curve  gives 
the  temperature  at  point  1  when  region  II  does  not  exist. 
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VI.  RESPONSE  OF  NUCLEI  TO  FIELD  FLUCTUATIONS 


A.  INTRODUCTION 

One  of  the  early  concepts  introduced  for  developing  a  y-Tay  laser  (Ref.  38)  was 
based  cm  die  idea  of  direct  emission  from  isomoic  nuclear  levels.  Such  levels,  because  of 
their  long  lifetime,  permitted  long  pumping  times  and  were  considered  to  be  good  storage 
levels.  Unfortunately,  because  of  the  narrow  natural  linewidths  associated  with  the  long 
lifetimes,  and  the  relatively  large  environmental  effects,  causing  inhomogeneous 
broadening,  resonance  conditions  for  transitions  involving  such  levels  are  difficult  to 
achieve.  Several  approaches  to  alleviate  the  problem  have  been  put  forth.  One  of  the 
approaches  proposes  to  reduce  the  line  broadening  by  external  fields  (Refs.  39  and  40).  It 
has  been  argued,  however,  that  in  order  to  narrow  the  inhotrH>geneously  broadened  line  by 
any  technique,  a  time  on  the  order  of  the  lifetinie  of  the  nuclear  state  would  have  to  be 
expended  (Ref.  27).  The  basis  of  the  argument  was  the  time  energy  uncertainty  principle 
(Refs.  41  and  42),  AEAt  >  h,  with  the  interpretation  that  AE  «  F,  the  natural  linewidth  of 
the  transition,  and  At  =x  the  lifetime.  It  was  stated  that  in  order  to  remove  inhomogeneous 
broadening  by  any  means  the  technique  used  would  have  to  be  applied  for  at  least  as  long  a 
time  as  t. 

We  have  argued  that  this  was  not  so  (Ref.  1,  Ch.  4)  and  that  in  fact  the  length  of 
time  the  line  narrowing  technique  would  have  to  be  applied  depended  (xi  the  strength  of  the 
interaction  causing  the  inhomogeneous  broadening.  The  appropriate  condition,  in  fact,  is 
that  the  time  for  narrowing  would  have  to  be  shorter  than  h/Heff  (Ref.  43,  p.  91 1)  where 
Heff  is  the  average  perturbation  causing  the  broadening.  The  time  energy  uncertainty  with 
AE  =  r  and  Ax  =  x  (as  interpreted  by  the  author  of  Ref.  27)  had  nothing  to  do  with  the 
phenomena  in  question.  The  time  required  to  remove  the  field  depends  on  the  strength  of 
the  field  and  not  on  the  ultimate  linewidth,  because  the  effect  of  the  external  fields  is  to 
cancel  out  the  inhomogeneous  field  by  shifting  the  resonance  lines  and  does  not  affect  the 
inherent  resolution  of  the  nuclear  energy  systenx 

In  this  paper  we  calculate  the  time  required  to  remove  or  set  up  a  hyperfine  pattern 
which  for  different  strengths  at  different  lattice  points  in  the  sample  produces 
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inhomogeneous  broadening.  The  calculation  is  based  on  the  Blume  theory  of  time- 
dependent  hypertine  interaction  and  was  originally  applied  to  relaxation  phenomena 
(Ref.  43). 


B.  THEORETICAL  BACKGROUND^ 


The  probability  of  emission  of  a  photon  with  energy  E  =  hco  and  wave  vectOT  k  by  a 
nucleus  making  a  transition  finm  state  e  to  g  is  given  by 


Peg(k.<0)  = 


(ho)+Eg-E/  +  (r/2)^ 


(35) 


where  is  the  interaction  Hamiltonian  between  the  nucleus  and  the  electromagnetic  field 
for  the  emission  of  the  photon,  le>  and  lg>  are  the  excited  and  ground  states  for  the 
nucleus  with  unperturbed  energies  Ee  and  Eg,  respectively,  as  shown  in  Hg.  14(a)  and  F  is 
the  total  linewidth  at  half  maximum.  For  tl^  rest  of  this  paper  we  assume  units  with  h  =  1. 


The  Hamiltonian  for  the  matter  system  including  all  the  nuclear,  atomic,  and  solid- 
state  terms  is  given  by  //,  so  diat 

H  le>  =  E,k> 


//  |g>  =  Eglg>  . 


To  obtain  the  experimentally  observed  lineshape  1(0))  we  have  to  average  over  all  the 
possible  initial  states  le>  and  sum  over  all  the  firud  states  |g>.  In  order  to  take  into  account 
the  dynamic  phenomena  due  to  time-dependent  Hamiltonians,  we  transform  to  the  time 
domain  to  get 


I«B) 


=  J  dt'  J  dt"  <  y^'V)  V^'^V)  > 


(36) 


where 


II 

1 

1 - - — 

a; 

1 

V  exp 

t 

-i  f//(t’)dt' 

1 

L  «  J 

V 

O 

(37) 


9  This  discussion  follows  the  analysis  of  M.  Blume  pfesenled  in  Ref.  43. 
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Figure  14(a).  Energy  levels  of  an  Fe^^  nucleus  without  a  magnetic  field  present 
(left)  and  with  a  magnetic  field  present,  causing  a  Zeeman  splitting  in 
the  levels  (right);  (b)  Absorption  spectra  of  Fe^^  corresponding 
to  the  level  splitting  shown  in  (a). 
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In  case  the  Hamiltonian  is  not  explicitly  time-dependent,  eq.  (37)  reduces  to 

V*(t)  =  e"‘ e-"‘  .  (38) 

and  if  in  the  experiment  only  photons  emitted  between  times  t  and  x  -i-  Ax  are  recorded,  the 
two  integrals  in  the  lineshape  eq.  (36)  should  be  taken  between  x  and  x  +  Ax  instead  of  0 
and  oo. 

In  the  case  of  a  standard  Mdssbauer  experiment  with  the  nuclei  unpolarized  and 
when  no  measurenoent  of  the  initial  state  is  made 

<  (t")  V"  (f)  >  =  <  (f  - 1")  >  , 

and  eq.  (36)  can  be  written  as 

-  r 

I«o)  =  -^  Re  J  dx  e  ^  <V  V  (x)  >  (39) 

O 

Let  us  consider  a  nucleus  in  a  time-dependent  magnetic  field  H(t)  =  h  f(t),  where  h 
is  die  maximum  value  of  the  magnetic  field  and  fit)  a  function  of  time  that  varies  between  0 
andl.  The  Hamiltonian  is  given  by 

//(t)  =  //o  +  gHl2hf(t)  .  (40) 

with  g  the  gyromagnetic  ratio,  p.  the  magnetic  moment  of  the  nucleus,  and  Iz  the  spin 
component  in  the  z  direction  which  is  also  the  field  direction.  Hq  is  the  ui^ierturbed  nuclear 
Hamiltonian  so  that,  Hq  ilimi>  =  Ei  llimi>  and  Hq  lloino>  =  Eo  llo™o>  with  Ii,  mi  and 
lomo  representing  the  excited  state  spin  and  z-component  (magnetic  quantum  number)  of 
the  nuclear  excited  and  ground  states,  respectively. 

Ccmsider  the  special  case  of  Fc57  with  Ii  =  3/2  and  Iq  =  1/2.  The  energy  levels  are 
shown  in  Fig.  14(a)  for  the  zero  field  case  on  the  left  and  for  the  nucleus  in  a  finite  field  H 
on  the  right.  Figure  14(b)  shows  the  corresponding  Mdssbauer  spectra  taken  in  the 
scattering  geometry  (total  scattoing  cross  section  measured). 

Suppose  the  magnetic  field  varies  with  time.  Blume  has  shown  that  in  this  case  the 
coirelatitxi  function  can  be  represented  as 


(41) 


C.  RESULTS 

Let  us  consider,  for  simplicity,  only  the  outer  lines  of  the  split  spectrum  and  the 
transition  (Iimi)  =  (3/2,  3/2)  to  (lonio)  =  (1/2, 1/2)  as  shown  by  the  solid  vertical  line  in 
Fig.  15. 

Equadtm  (41)  dien  reduces  to 

t 

-i6  f  f(t')dt‘ 

<V  V  (t)>  =  e  (5(3/2, 1/2, 1/2, 1/2) e  “  ,  (42) 

with 

5  =  (g,mj  -  goHio)  nh  and  *  Ej  -  . 

and  Q  is  die  time  independent  transition  matrix  dement 

If  we  further  assume  that  the  field  varies  exponentially  with  t  with  time  constant  y 
and  shift  the  energy  reference  by  ©lo,  then 

f(t)  =  c’^*  . 

o 

and  the  part  of  the  spectrum  representing  the  (Iimi)  =  (3/2,  3/2)  ->  (lomo)  =  (1/2, 1/2) 
transition,  leaving  out  the  factor  Q,  is  given  by 


I(©)  =|r  Re  J  dt  exp|^-i  tor  -  (Tnyt  -y-  (l  -  c  (43) 

O 

This  can  be  integrated  numerically  to  obtain  the  lineshape  for  arbitrary  values  of  the 
parameters.  Fen’  two  limiting  cases  the  solutions  can  be  expressed  in  closed  form. 
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Figurw  IS.  Absorption  cross  ssetion  of  Fs^^  showing  only  the  outer  lines, 
ss  indiested  in  Fig.  14(b)  when  r  s  1.0  snd  5  s  1.0  snd  y  s  0.1  in  (a) 
end  Y  s  KH  in  (b)  snd  y  »  5  x  10^  In  (e) 


Case  (a)  when  F-  5  and  y»  F  we  get 


1(0)) “I?  Re  ^  y -itoj  y  +7-1(0^ 


^  -iet)j  y  +Y-itDj  y  +2y-u 


8  (y+  F)  (d 


fl^  2 

2  V  +  0)^ 

(4 


fr"  ]\(r  f  * 

+0)2^ 


5=  w(^r+ir].£ 

.(t  .(f 


Clase  (b)  when  F  »  Yand  F •  8  we  get 


=  |.Re 


where 


iRc-1 


^-i((a  +  8)  ^  (-1)”*  (2m)t 

(to  +  5)^  +  ■^—  ®“® 

4 
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The  first  term  in  the  series  is  1  so  that  in  die  limit  F  »  Ythe  line  shape  approaches  a 
Lprentzian  cratered  at  A(o  and  with  a  full  width  at  half  maximuin  of  F. 

Figure  15(a)  shows  a  split  spectrum  for  values  of  y  <  F,  and  5  =  F  (y  =  10“^ 
F  =  1, 5  «  1).  Figure  15(b)  shows  a  collapsed  spectrum  for  y  =  F  =  5  and  Fig.  15(c) 


shows  an  enhanced  collapse  with  a  narrower  line  spectrum  for  y  =  Sr  and  5  =  F. 
Figure  16(a)  shows  a  completely  split  spectrum  with  y/T  =  10~^  and  S/F  =  10,  while 
Fig.  16(b)  shows  a  completely  narrowed  single  line  with  y/F  =  lO^  and  8/F  =  100. 
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Figure  16.  Absorption  cross  section  of  showing  oniy  the  outer  lines  as 
indicated  In  Fig.  15(b)  for  the  case  when  y  &  10^,  r  s  i.o,  5  s  10  in  (a), 
and  Y  s  IQS,  r  s  1.  8  B  100  in  (b). 


These  results  show  that  the  nuclear  resonant  lineshape  will  show  a  full  split 
spectrum  as  in  Fig.  16(a)  if  the  change  of  field  is  slow  (y  «  O  compared  to  the  decay  time 
of  the  nucleus.  In  this  situation  most  nuclei  will  be  exposed  to  the  field  during  their  lifetime 
and  while  they  decay.  If  the  field  decay  is  fast  (y  »  F)  the  spectrum 
will  show  a  collapsed  lineshape  similar  to  the  one  shown  in  Fig.  14(b)  on  the  left.  This 
implies  that  it  does  not  require  time  on  the  order  of  a  lifetime  of  the  nuclear  state  to 
effectively  remove  the  inhomogeneous  broadening. 

D.  CONCLUSIONS 

Special  cases  for  low  and  high  relaxation  y  «  F  and  y  »  F,  respectively,  with 
F  »  5,  indicate  that  the  natural  lifetime  of  the  nuclear  state  F  sets  the  lower  limit  on  the 
width  of  die  spectra  obtained  for  these  cases  but  does  not  for  fast  relaxation  prevent  a  good 
definition  of  the  spectrum.  Thus  in  a  very  shot  time  during  which  the  hyperfine  interactitm 
changes  value  the  nuclear  resonance  spectrum  follows  this  change  quickly.  In  a  time  much 
less  than  the  natural  lifetime  the  spectrum  changes  fiom  being  rq}resentative  of  one 
hyperfine  interaction  to  being  representative  of  another.  The  nuclear  system  quickly 
respcmds  to  the  change  of  fields  in  its  environment  Thus,  if  inhomogeneous  broadening  is 
represented  by  a  sum  of  spectra  representative  of  different  fields  at  different  nuclear  sites,  if 
these  fields  can  be  changed  in  a  time  At  much  shorter  than  the  natural  lifetime,  the  resulting 
spectra  taken  over  a  lifetime  of  the  nuclear  state  will  represoit  the  final  hyperfine  interaction 
spectrum. 

This  result  is  important  to  the  y-ray  laser  feasibility  study  because  it  shows  that 
inhomogeneous  broadening  can  be  reduced  quickly  by  applying  sufficienuy  strong  and 
appropriate  external  fields.  Inhomogeneous  broadening  of  several  (»ders  of  magnitude 
higher  than  the  natural  lifetime  of  the  isomer  can  be  removed  by  fields  in  a  time  shorter  than 
the  natural  lifetime  of  the  isomer.  Such  a  fast  mechanism  could  be  used  to  trigger  the  lasing 
or  supetfluorescent  emission  by  providing  a  strong  resonance  overlap  of  narrow  resonance 
lines  quickly. 
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APPENDIX  A 

DERIVATION  OF  THE  HAAKE-REIBOLD  EQUATIONS 


The  purpose  of  this  appendix  is  to  justify  the  system  of  eqs.  (1-3)  in  Chapter  I  for 
the  excited  population  number  operators  Ni,  the  polarization  operators  Rij^,  and  the  electric 
field  operators  These  equations  are  given  in  Ref.  A-1  without  derivation. 

Reference  A-2  derives  similar  equations  for  a  two-level  atmnic  system,  but  does  not 
indicate  how  to  generalize  the  theory  to  multi-level  systems.  This  appendix  will  attempt  to 
fill  in  the  steps  needed  to  extend  the  derivatitx)  of  Ref.  A-2  to  die  present  case  of  interest. 

For  this  purpose  it  will  be  necessary  to  use  generalized  commutation  relations  given 
by  Ref.  A-3  for  operators  Skj*  that  effect  a  transition  from  state  j  to  state  k  in  a  multi-level 
system.  The  commutation  teladtxis  are 

[Sjk*  Smnj  =  Sjn^km  ‘  ^ink\j  * 


Reference  A-3  neither  derives  nor  gives  a  reference  for  (A-1),  which  does,  in  fact, 
reduce  to  the  familiar  two-level  commutation  relatitxis  (cf..  Ref.  A-4,  which  derives  them) 

[r*.R-].  R3.[R3-''‘]=±R*  (A-2) 

for  the  polarization  R^  and  inversion  R3.  The  identification  between  the  operators  labeled 
S  in  (A-1)  and  the  operators  labeled  R  in  (A-2)  required  for  this  reducticxi  is 

Sjo  =  R  ♦  Sqj  =  R  ,  Sjj  -  Sqq  =  R^  .  (A-3) 


*  To  avoid  confusion  with  the  electric  field  the  letter  S  is  used  here  to  designate  the  transition  operators 
instead  of  E,  which  Ref.  A-3  uses  for  that  purpose. 
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The  identification  eq.  (A-3)  for  the  two-level  case  suggests  the  appropriate 
corresponding  identification  for  the  operates  pertinent  to  the  multi-level  case  leading  to  the 
eqs.  (1-3)  in  Chapter  I.  Thus,  generalizing  eq.  (A-3),  it  will  be  assumed  that 

Sy  =  Nj,i  =  0, 1,2,3  , 

Sii  .Sj-Ri:.i  =  2.3,j  =  1.2,  i>j  .  (A4) 

Following  Ref.  A-2,  the  Hamiltonian  for  the  non-interacting  atoms  is  given  by  a 
sum,  over  all  amms  in  the  system,  of  terms 

=  h  j+  0)32[N3-N2]+  j"*"  “43(^4‘^3  ]■*■  “lo(^r^o)  *  (^-5) 

The  interaction  Hamiltonian  is  given  by  a  sum,  over  all  atoms  in  the  system,  of  terms 

H„  =  ih  L  E,,  [R,:  -  *3^  )•  gjEjj  [R,!  -  I  .  <A-6) 

The  Heisenberg  equation  of  motion  for  any  operator  0  of  interest  in  terms  of  the 
Hamiltonian  H  will  have  die  form 

dO/dt  =  (i/h)  [H,  0]  .  (A-7) 

However,  tite  eqs.  (1-3)  of  Chapter  I  cannot  be  derived  entirely  from  a  relation  of  the  form 
eq.  (A-7);  the  loss  and  source  terms  are  ad  hoc  additions  to  equations  that  are  obtainable 
from  a  Hamiltonian. 

The  Hamiltoiuan  to  be  used  to  obtain  the  Heisenberg  equations  before  the  addition 
of  the  od  hoc  terms  consists  of  the  standard  free  photon  field  contribution  plus  the  sum  over 
all  atoms  of  the  single  atom  expressions  given  by  eq.  (A-S)  and  eq.  (A-6).  On  substituting 
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the  Ni  and  the  Rij*  into  the  resulting  equation  of  the  form  eq.  (A-7)  and  then  using 
eqs.  (A-1,  A-3,  A-4),  the  following  equations  are  obtained: 

dNj/dt  =  gjE3j[R31  +  *^2]  . 
dN2/dt  =  g2E32[^3l'*‘^2]  . 

dN3/dt  =  -[g,E3j(R^l  +  R^l]+g2E32(l^32+^2]]  . 

dR*/dt  =  i  (±  “31  +  “32  ±  ^12  +  -  “lo]  1^1  +  +  g2E32R*, . 

dR^j/dt  =  i  “32  *  ^“31  +  “21)  +  gj  (^3  "  E3J  -  g2E32R21  ’ 

dR^2/<lt  =  i  “31  ^  20)32  +  0)2j]  (N3  -  N2]  _  g,E3jR2i  . 

Equations  (A-8)  are  at  the  stage  of  conesponding  equations  for  the  two-level  atomic 
system  in  Ref.  A-2  which  Ref.  A-2  then  modifies  by  separating  the  electric  field  and 
polarization  operators  into  their  positive  and  negative  frequency  parts,  neglecting  rapidly 
oscillating  terms  in  the  process.  Applying  the  same  process  to  (A-8)  results  in  the 
eqs.  (1-2)  of  Chapter  I  minus  the  source  and  loss  terms.  Equation  (la)  is  not  included 
because  the  right  hand  side  consists  only  of  a  source  and  a  loss  term. 

Equations  (3a,  3b)  in  Chapter  I  can  be  derived  from  Maxwell's  equations  with  a 
polarization  source  exactly  as  in  Ref.  A-2  by  assuming  that  the  radiation  due  to  a  particular 
transition  depends  only  on  the  polarization  associated  with  that  transition.  The  results 
appear  to  differ  somewhat  in  form  from  the  corresponding  results  given  in  Ref.  A-2. 
However,  it  can  be  seen  by  inspection  that  the  differences  are  due  to  the  use  of  a  retarded 
time  coordinate  and  neglect  of  spectral  broadening  in  Chapter  I  as  in  Ref.  1. 

Actually,  two  discrepancies  between  the  equations  derived  in  the  manner 
summarized  in  this  appendix  and  eqs.  (1-3)  wiU  exist  as  eq.  (A-8)  now  stands.  These  are 
the  presence  of  the  coupling  constant  gi  in  aU  six  equations  and  the  presence  of  a  term 
involving  the  oii  in  each  of  the  last  three  equations.  The  first  discrepancy  can  be  accounted 
for  as  due  to  a  normalizatitHi  of  the  electric  field  and  the  second  as  due  to  irrelevant  phase 
factors  in  die  polarization  operators. 
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COMPUTER  CODE  FOR 
TRIPLE  INTEGRAL  CALCULATION 


Bl.  BACKGROUND 

In  this  appendix  we  describe  the  algorithm  used  to  compute  the  superfluorsecent 
pulse  from  the  closed  form  (triple  integral)  expression  derived  in  Chapter  1.  The  physics  of 
the  problem  is  modeled  by  the  set  of  differential  equations  with  stochastic  source  terms 
given  in  chapto’  I.  section  2. 

The  algorithm  makes  use  of  the  following  expressions  derived  in  chapter  I: 

Equation  (1): 


Equation  (2): 


i  r  'f’  -T 

k(t,t’)  =  2  j  N3(t")  dt"  +  No  ^  y  — J  -  No  (t  -  f) 


Equation  (2*): 


k(U')  =  ko(t)-ko(t') 


C  27N0  cil  c21  lei: 
i.e.,  ko(t)--=r—  1  r  T  J'^ol  Y  * 


and  ko(t)  =  2  J  NjCt")  dt"  -  Nq  * 


Equation  (3): 


B2.  STEP  BY  STEP  DESCRIPTION  OF  ALGORITHM 

Fig.  B1  shows  the  step  by  step  iteration  process  used  in  the  algorithm.  Imtially  the 
linear  approximation,  function  I(t),  is  used  to  calculate  two  starting  pointsby  iteration;  then 
from  these  a  two  point,  linear,  extrapolation  procedure  is  set  up  and  used  to  project  an 
initial  point  for  the  next  set  of  iterations.  The  steps  are  indicated  and  numbered  on  the  left, 
while  the  equations  from  which  the  equations  in  the  codes  are  derviced  are  indicated  on  the 
right  in  paranthesis. 

Step  1  Calculate  f(t)  on  the  range  t »  0  to  t  =  tn  from  (1)  and  (5) 
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Figure  Bl.  Schematic  showing  numerical  procedure  for  calculating  SF  pulse. 

I  (t)  gives  the  linear  solution  and  l(t)  the  final  solution,  Including  the  non-llnear  terms. 


Step  2 


Use  t(t)  to  calculate  Ii(t)  and  Ii(ti)  from  equations  (2),  (3)  and  (1)  (with  I(t) 
substituted  for  t(t)): 


N3(to)  = 


YNq  f  -Yt 


e-^)-P  t(t’)  df 


(3) 


where  P  =  i  +  |il 


k(U‘)*2 


t 

jN3(t")dt"  +  No 
** 

e  -  e 

Y 

t' 


No(t-f) 


(2) 


Ii(to)  =  ^c-^‘ IJ  ^[2^k(U’)  x’]e-'“’''^’^^*'dx'dt’ 

Step  3  Use  points  Ii(to)  and  Ii(ti)  to  calculate  slope  Soi  of  line  I  oi,  and  then  use  Sqi 
and  At  =  t2-ti  to  calculate  I2(t2)  which  is  initial  estimate  for  I(t2). 

Step  4  Use  eqs  (1),  (2)  and  (3)  to  calculate  first  iteration  estimate  for  I(t2)  in  the 
following  way:  use  I(t2)  for  I(t)  in  eq.  3  and  then  the  resulting  N3  in  eq.  2  to 
get  first  iteration  result  for  I(t). 

Second  iteration  -  use  I(t)  in  step  3  and  4  to  calculate  the  new  I(t).  Then,  with 
these  two  point  replacing  Ii(to)  and  Ii(ti),  repeat  step  2  and  continue  this 
process  until  convergence  is  satisfactory. 


RELATIVE  INTENSITY 


>— eg  t.  - 1— _ I _ I _ I _ I _ I 

0.000  1.238  2.478  3.713  4.950  6.188  7.425  8.663  9.900 


TIME 


2-1-904 


Figure  B2.  Sample  output  for  N  s  10^  and  y  s  50,  r  s  10*^  and  At  s  T  a  ti4.i  - 1|  s 
0.1,  showing  the  linear  solution  r(t)  and  the  final  reauh  l(t).  Time  Is  In  units  of 


B- 


RELATIVE  INTENSITY 


0.064 

0.056 

0.048 

0.040 

0.032 

0.024 

0.016 

0.008 

0.000 

0. 

Figure  B3.  Another  eample  run  with  N  s  10'*°  showing  the  finai  resuit  after  an 

average  of  three  iterations  per  time  point.  # 
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DESCRIPTION/DATA  AND  VARIABLE  NAMES 


Following  is  a  description  of  the  data  written  to  the  various  unit  numbers  in  the 
FORTRAN  program.  In  the  command  procedure  these  unit  numbers  are  assigned  to 
ASCn  files  which  can  be  edited  or  printed. 


I 


) 


IM  Description  of  the  Data 

12  T,  iteration  #1  to  3;  N3,  result  of  eq.  1  w/o  the  dt’  integral; 

,  result  of  eq.  1 

14  1st  30  entries  in  the  k(t,t')  array 

15  T;  eq.  I  before  integral;  integral  in  eq.  1,  eq.  1,  counter 

17  T,  eq.  1  before  double  integ^,  Na;  Nb;  N3,  result  of  eq.  1  w/o  the 

dt  integral,  slope  between  Ais  poin  t  and  previous  point,  slope  value 

Variable  Names 


Text 

Computer  Code 

g 

G 

Y 

GAMSML 

(input) 

N 

N 

(input) 

r 

GAMLRG 

(input) 

MU 

(input) 

No 

NO 

(input) 

2* 

Besl 

-1 

1 

->< 

GLMGS 

Y-r 

GSMGL 

/ 

LBAR 
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GSTGSDN 


N 


dx  increment 


location  in 


I 

J 


minimum  t 


t 

in  J(usually  0) 


thet 


inj 


increment  in 


J 


O 


DELTAX  (input) 

TCHANGE 

T 

TMIN 

TMAX 

TINCR 


The  subroutine  name  and  equation  number  in  the  text  are  correlated  in  the  following  way:  • 

Subroutine  Eq.  1  is  equation  number  1. 


B-10 


LISTING  OF  PROGRAM 


PRCX3RAM  BESAPRJ2 
Written  by  Karen  K.  Garcia 
A]^  21. 1989 

This  program  uses  a  Bessel  function 
from  the  IMSL  routines 
PARAMETER  (NS1ZE=1000JCS1ZE=1000) 

DIMENSION  IPK1(300) 

REAL  X(NSIZE),  EQ16ARR(NSI2E),  NO.  MU.  G.  LBAR. 

N.  N3.  KTTP(KSIZEJCSIZE).  NA.  NB.  NAPARTl 
CX)MMON  /BLOCKl/  NO.  MU.  GAMLRG.  GAMSML.  KTTP 
COMMON  /BLOCK2/  YINPT.  G.  N3.  Y(NSIZE).  Y2(NSI2E).  NA.  NB 
COMMON  /BLOCKS/ NAPARTl.  BTA,  SL0PEVAL(NS1ZE).  TCHANGE 
COMMON  /BLOCK4/  EQ20T1 
G  =  1.0 

WRITE(6.*)'Enter  min.  max  and  incr  of  time* 
READ(5.*)TMIN.TMAX.TINCR.tmax2 

WRrrE(6.*y  ist  tmax=  ‘.tmax 
WRnE(6.*)*  1st  TMAX2*  •.TMAX2 

WRITE(6.*)'Enter  delta  x  value* 

READ(5.*)DELTAX 

Read  in  the  values  fw  the  input  data 
WRITE(6.*)TVpe  value  for  No* 

READ(5.*)NO 

WRITE(6.*)Type  the  value  for  MIT 
READ(5.*)MU 

WRlTE(6.*)Type  the  value  for  large  gamma* 

READ(5.*)GAMLRG 
WRrTE(6.*)Type  the  value  for  N* 

READ(5.*)N 

Check  the  initial  value  for  Ibar 
Ibar  B  0.1 
1CNT=  1 


WRITE(6.*)Type  the  value  fca  small  gamma* 
READ(S.*)GAMSML 

WRrrE(6.*)Type  value  for  the  start  of  the  slrqre’ 
READ(5.*)TCHANGE 

GLMGS  B  GAMLRG  •  GAMSML 

GSMGL  B  GAMSML  •  GAMLRG 

GSTGSDN  B  (G**2  *  GAMSML)  /  N 

NAPARTl  B  (GAMSML  *  NO)  /  (GAMLRG  -  GAMSML) 

BTA  «  (1.0  /  G)  +  (2.0  •  0.5*MU  •  LBAR) 

EQ20T1  =  (2.0  •  GAMSML  *  NO)  /  GLMGS 


YMINB9.0e20 

YMAXb-9.0E20 

Ib  1 

NCALLbI 

OPEN(UNIT»15J>ISPOSE»*KEEF.STATUSb’NEW) 

WRnE(15.15) 

FORMAT(2x.T.4X.'EQ16TERM  l’.2X. 

'  INTEG’.2X.'EQ16&20*3X.TCN1-) 

Check  to  see  if  tmin  B  tmax 


on  on 


IF(TMAX-TMINX.T.TINCR)GOTO  25 
ITCRsO 
T»TMIN 

DO  WHIIJE(T1£.TMAX) 

CALL  EQ1(TSUM,TINCRJ4CAUL,TJTER. 

ICNT,GLMGSJ>ELTAX) 

Save  the  X  values 
X(I)  =  T 

Compute  die  equation  labelled  16 
EQ16T1  =  GSTGSDN  *  EXP(-GAMLRG  •  T) 

EQ16ARR(I)  =  EQ16T1 
Y(I)  *  EQ16T1  *  TSUM 
YMIN  =  MINl(YMIN.Ya)) 

YMAX  =  MAX1(YMAX.Y(I)) 

WRnE(15^)X(I)JEQlCTl.TSUM,Y(04 
20  FORMAT(1XJF6.2^X^^X£10.43X310.43X44) 

I»I+  1 

C  Check  against  the  airay  size 
IFa.GTJ4SIZE)THEN 

WRnE(6.*)TNCREASE  ARRAY  DIMENSION.  (NSIZE)’ 

CALL  EXIT 
END  IF 
T  *  T  +  TTNCR 
END  DO 

25  CONTINUE 
CLOSE(UNrr»15J)ISPOSE=’KEEF) 
OPEN(UNrr«14J)ISPOSE»XEEP^ATUS»’NEWO 

C  write  out  the  KTTP  array 

wrilB(14,*)  'columns  1  •  10' 
write(14^7) 

27  foniiat(lx,10x.'0',10x.'l',10x,'2',10ji,3',10x,'etc.') 

DOn«130 

WRnE(14^ii.(KTrP(njJ)JJ=l,10) 

END  DO 
wrile(14.*)" 

write(14  *)  'columns  11  •  20' 
wrile(14*)" 

DOn»130 

WRITE(14^.(KTTP(njJ)JJ=11^0) 

END  DO 
wrilc(14,*)' ' 

writB(14  *)  'columns  21  -  30' 
wrilc(14*)" 

DOnsl30 

WRITE(14.2Qii.(KTIP(njJ)JJ=2l30) 

END  DO 

CLOSE(UNIT«14J>ISPOSE«'KEEP') 

26  FORMAT(i4.10(lX310.4)) 

C  END  OF  FIRST  EQUATIONS  SET.  START  OF  SECOND  EQUATIONS  SET 

^**********«***«********»****4i****«i***4t**»****«t*****4i**»»**»**»4i****»»4 

c 

I«I.l 

LASll-I 

NCALL>2 
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c 

C  Write  the  headers  for  the  second  equation  set 

OPEN(UNIT=12J)ISPOSE='KEEP,STATUS=’NEWl 

OPEN(UNrr»17J)ISPOSE=*KEEF,STATUS=’NEW*) 

WRnE(12.30) 

30  FORMAT(5X.T,5X.’J'.8X,'N3M5X.THAT(T)*. 

*  12X,TPRIME(t)') 

WR1TE(17.199) 

199  F0RMAT(4X,T;2X,'EQ1  term  r.4X.'NA*.7X,W. 

*  7X.'N3’.6X.'Ihat(t)'.4X;Slope  value’.4X. 

*  l(t)  Nth’^X.’SLOPE') 

C  Compute  die  second  set  ot  equatkxis 

ICNTsO 
C 

C  6/6/89KKGCopy  the  Diat  array  to  the  new  I  array 

C  will  be  used  fw  the  iterations 

DOIFlLL=  l^ASTI 
Y2(IFILL)  =  Y(IFILL) 

END  DO 
C 

DOJ=  1,LASTI 

C  Restore  the  x  value  from  the  above  computation 
T  =  X(J) 

C 

C  Iterate  over  the  re-substitution  of  die  I(t)  term. 

DOK«U 

C 

C  Equation  1  parameters 

C  TSUM  -  output 

C  TTNCR.  NCALL.  T  -  input 

CALL  EQlCrSUM.TINCRJfCALL.Tj:, 

*  ICNT.GLMGSJ)ELTAX) 

C 

YINPT  =  EQIOARRO)  *  TSUM 
Y2(J)  =  YINPT 

WRnE(1235)TJLN3,Y(I).YINPT 
35  PORMAT(1X46.2.1X.I4.3(1X£20.13)) 

END  DO 
IF(j.GT.l)THEN 
C 

C  Compute  the  sk^ 

SLOP*(Y2(j)-Y2(j-l)) 

*  /(X(D-X(i-1)) 

WRITE(12,*)'SLOPE  BETWEEN  ’j-l, 

*  •  AND  '  j.’  IS  ’.SLOP 

write(12,*)’the  values  are  y2(j  j-l).x(j  j-l)',y20'), 

*  y20-l).xa).x(j-l) 

END  IF 

IF(T.GE.TCHANGE)THEN 
SLOPEVAL(J+l)  *  Y2(J)  +  (SLOP  *  TTNCR) 

Y2(J+1)  =  Y2(J)  +  (SLOP  •  TTNCR) 

END  IF 

C  write  data  to  unit  17 

IF(J.EQ.I)THEN 

WRnE(I7.200)T£Ql6ARR(J)JfAJfBJO.Ya). 

*  SLOPEVAL(J+I).Y2(I) 
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ELSE 

WRITE(17^00)T£Q16ARR(J);^AJ4B;^3,Ya). 
•  SLOPEVAL(J+l).Y2(J).SLOP 

END  IF 

200  FORMAT(lXJJ6^^(2X£8i).4(3X£10.4)) 

END  DO 

aX)SE(UNrr*12DISPOSE*'KEEP) 

CLOSE(UNrr=17J)ISPOS&='KEEP) 

CKG  END  DO 

END 


ALL  SUBROUTINES  BELOW  HERE 


EQUATION  1 


SUBROUTINE  EQ1(TSUM.TINCRJ^CALL,TJTER. 

*  ICNT.GLMGSJ)ELTAX) 

Equation  1  i)aiaineters 
TSUM  -  output 
TINCR,  NCALL,  T  •  input 

PARAMETER  (NSIZEsl000JCSIZE=1000) 

REAL  NO.  MU,  N.  N3.  KTTP(KSI2EJCSIZE)JJAJ4B.G.NAPART1, 

*  TEMPEQ16T2,  EQ1CT2 

DOUBLE  PREQSION  EQ16T3,  EQ16T3N2.  DEQ20.  DXPRIME 
CX)MMON  /BLOCKl/  NO,  MU.  GAMLRG.  GAMSML,  KTTP 
COMMON  /B1JOCK2/  YINPT.  G.  N3.  Y(NSIZE),  Y2(NSIZE). 

*  NAJffi 

COMMON  /BLOCKS/  NAPARTl,  BTA,  SLOPEVAL(NSIZE).  TCHANGE 
COMMON  /BLOCK4/  EQ20T1 
C 

C  VARIABLE  DEPINmON 

C  NCALL  •  Number  of  die  equation  set  that 

C  this  subioutine  is  called  for. 

C  ix.  1  -  forGrst  curve 

C  2 -for  second  curve 

C 

TrSUM»0.0 
TSUM  «  0.0 
TNEWSUM  s  0.0 
C  Check 

IFCriT.Tincr)GOTO  20 
TPRIME  »  0.0  • 

TLAST»T-TWCR 

jOflgBO 

DO  WHILEGTRIMEiE.TLAS'D 
C 

IF(NCALLJEQ.1)1HEN 

IBSLFLG«0 

CALL  EQ5(EQ20,TPRIME.T,TlNCR,GLMGS,tppt) 
IF(EQ20JLT.0.0)THEN 

WRnE(6.*yJ0  USED  Isi  eq.  EQ20b  •EQ20.*  tprime*  '^prime 
mSLFLG-l 
END  IF 

ELSE  IF  (NCALLEQ.2)THEN 
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c 


mSLFLGaO 

CALL  EQ2CnNCR£(320.TPRIME.TJT1ER4BSLFLG) 

END  IF 

XPRIME»0.0 
TXSUM  =  0.0 
KCNT=  1 

XLAST=  1.0-XPRIME 
EQ1CT3P2  =  (GLMGS  •  TPRIME) 

DO  WHILE(XPRIMEXE.1.0) 

Compute  the  equation  labelled  16  integration  portion 

XPPT  =  XPRIME  +  DELTAX 
IF(IBSLFLG£Q.O)THEN 
DEQ20sEQ20 
DXPRIME  a  XFRIME 

EQ16T2  =  DBSI0(2.0D0  •  DSQRT(DEQ20  *  DXKUME)) 
ELSE  IF(IBSLFLG£Q.1)THEN 
IF(J0FLG£Q.0)THEN 
WRnE(6  *)*J0  BESSEL  FUNCHON  USED’ 
WRITE(6*)*XPRIME,THUME'^RIME,TraiME 
J0FLG=  1 
END  IF 

7/3/89  Put  Itv  Kay's  JO  bessel  function  in  here 
BESL  =  (2.0  •  SQRT(EQ20  *  XPRIME)) 

POWER  SERIES 
IF(BESL.LT.10.0)THEN 
DONBES»1.18 

FITOP  =  (-((BESU2.0)**2)**NBES) 

Compute  N  factorial 
rFACr=1.0 
DOL  =  NBES.l,-l 
rFACT  =  iFACT  *  float(L) 

END  DO 

FlBOT  =  rFACT**2 
EQ16T2  *  (FITOP  /  FIBOT)  +  1.0 
END  DO 

ASYMPTOTIC 

ELSE  IF(BESL.GE.10.0)THEN 
PI  =  3.141592654 

F2TRM1  =  SQRT(2.0  /  (PI  *  BESL)) 

F2TRM2  =  COS(BESL  -  (PI  /  4.0)) 

EQ16T2  =  F2TRM1  *  F2TRM2 
ENDIF 
END  IF 

EQ16T2  =  EQI6T2**2 
EQ16T3N2  =  0.0D0 
EQ16T3  =  0.0D0 

EQ16T3P1  =  -2.0  •  0.5*MU  *  XPRIME 
EQ16T3  =  EXP(EQ1CT3P1  +  EQ16T3P2) 

EQ16T3N2  «  EXP(EQ16T3P1  +  GLMGS  •  TPPT) 

EQ16T3  » (EQ16T3  +  EQ16r3N2)  /  2.0D0 
IF(EQlOT3iT.1.0E-10)THEN 
EQ16T3  =  l.OE-10 
ENDIF 
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EQ16  =  EQ16T2*EQ1CT3 
IF(IXCNT£Q.1)THEN 
TXSUM  =  EQ16  *  DELTAX 
PREVTX  =  TXSUM 
IXCNT*IXCNT+1 
TXSUM  =  0.0 
ELSE 

TXSUM  =  TXSUM  +  ((EQ16*DELTAX)  +  PREVTX)/2.0 
PREVTX  =  EQ16  *  DELTAX 
END  IF 

XPRIME  =  XPRIME  +  DELTAX 
END  DO 
mSLFLGsO 
IFacnLEQ.DTHEN 
TTSUM  =  TXSUM  •  TTNCR 
PREVTT  =  TrSUM 
ICNT  =  ICNT+  1 
TTSUM  =  0.0 
ELSE 

TTSUM  =  TTSUM  +  ((TXSUM*TINCR)  +  PREVTT)/2.0 
PREVTT  =  TXSUM  •  TTNCR 
END  IF 

TPRIME  =  TPRIME  +  TTNCR 
END  DO 
TSUM  =  TTSUM 
20  CONTINUE 
C 

RETURN 

END 


EQUATION  2 


SUBROUTINE  EQ2(TINCR.EQ20,TPRIME.TJIER4BSLFLG) 
PARAMETER  (NSIZE=1000JCSIZE=1000) 

REAL  NO.  MU,  N,  N3,  N3SUMSAV,  N3SUMPREV,  N3SUMCUM, 

*  KTIPCKSEEJCSIZE),  KOT.  KOTPR,  NA,  NB  JC0TT2, 

*  KOTTPn,  NAPARTl 

COMMON  /BLOCK!/  NO,  MU.  GAMLRG,  GAMSML.  knp 
COMMON  /BLOCK2/  YINPT,  G,  N3.  Y(NSIZE).  Y2(NSIZE), 

*  NAJffi 

COMMON  /BLOCK3/  NAPARTl,  BTA,  SLOPEVAL(NSIZE).  TCHANGE 
C  Compute  equation  2 
XN3  =  0.0 
xn3sum  =  0.0 
IFCriT.TINCR)GOTO  20 
TVAL  =  0.0 

K0TT2  =  EXP{^AMSML  •  TT/GAMSML 
DO  WHILE(TVALiE.T) 

IF((TVAL-TINCR)iT.TINCR)THEN 
N3  =  0.0 
XN3  =  0.0 
PI-SF. 

TDPRIME  =  TVAL  -  TTNCR 
CALL  EQ3(TTNCR,TDPRIMEjrER) 

XN3«N3 
END  IF 

TDPRIME  »  TVAL 


B-17 


noon 


CALL  EQ3(TINCR.TDPRIME4TER) 

XN3  =  (XN3  +  N3)  /2.0  •  TTNCR 
xn3suin  =  xn3suin  +  xn3 
TVAL  =  TVAL  +  TTNCR 
END  DO 

20  CONTINUE 
xn3  =  xnTsum 

KOT  =  2.0  *  XN3  -  NO  *  (K0TT2  +  T) 

C 

XN3  =  0.0 
xnTsuin  =0.0 

IF(TPRIMEiT.TTNCR)GOTO  25 
TVAL  =  0.0 

K0TTPT2  =  EXP(-GAMSML  *  Tprime)/GAMSML 
DO  WHILE(TVALJJE.TPRIME) 
IF((TVAL-TTNCR)XT.TTNCR)THEN 
N3  =  0.0 
XN3  =  0.0 

PT 

TDPRIME  =  TVAL  -  TTNCR 
CALL  EQ3(TTNCR,TDPRIME.ITER) 

XN3  =  N3 
END  IF 

TDPRIME  =  TVAL 

CALL  EQ3CnNCR.TDPRIME4TER) 

XN3  =  (XN3  +  N3)/2.0  *  TTNCR 
xnSsum  =  xn3sum  +  xn3 
TVAL  =  TVAL  +  TTNCR 
END  DO 

25  CONTINUE 
xn3  =  xnTsum 

KOTPR  =  2.0  *  XN3  -  NO  *  (K0mPT2  +  Tprimc) 
EQ20  =  K0T-K0TPR 


C  5/30/89  KKG  Add  code  to  set  flag  to  use  JO  bessel 
c  function  if  eq20  less  than  0.0 
IF(EQ20XT.0.0)THEN 
IBSLFLG=1 
EQ20  =  ABS(EQ20) 

END  IF 
RETURN 
END 


EQUATIW3 


SUBROUTINE  EQ3(TTNCR,TDPRIME  JTER) 
PARAMETER(NSIZE=10004CSIZE=1000) 

REAL  NO,  MU,  N,  NBSUM,  N3,  NA,  NB, 

*  NAPARTl,  KTTP(KS1ZEJCSIZE),  NAPART2,  lOFT, 

*  NBSUMCUM,NBSUMSAV,NBSUMPREV,NB1,NB2 
COMMON  /BLOCKl/  NO,  MU,  GAMLRG,  GAMSML,  kttp 
COMMON  /BLOCK2/  YINPT,  G,  N3,  Y(NSIZE),  Y2(NSIZE), 

*  NA,NB 

COMMON  /BLOCK3/ NAPARTl,  BTA,  SLOPEVAL(NSlZE),  TCHANGE 
NA  =  0.0 
NB  =  0.0 

C  Compute  equation  3 
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NAPART2  *  EXP(-GAMSML  *  TDPRIME)-EXP(-GAMLRG  •  TDPRIME) 
NA  =  NAPARTl  •  NAPART2 
NBSUM  =  0.0 

IF(TDPRIMEl,T.TINCR)GOTO  20 
ICNT=  1 
TVAL  =  0.0 

DO  WHILE(TVALiE.TDPRIME) 

NBl  =  EXP(GAMLRG  *  (TDPRIME  -  Tval))  •  Y2aCNT) 
nbsiun  =  nbsum  +  (nbl'*'tincr) 

ICNT  =  ICNT+1 
TVAL  =  TVAL  +  TINC31 
END  DO 

NB  =  BTA*  NBSUM 
20  CONTINUE 
N3  =  NA  -  NB 
RETURN 
END 
C 

Q****m***0m*m***»**»***m********0****m0»**m*mm***^**m****m******** 

C  EQUATIONS 

Q000*0000*0m*0m0*000*0»*0*0m00**m00*0m00mm*mm**m00***»*»»*****0*** 

c 

SUBROUTINE  EQ5(EQ20.TPRIME.T.TINCR.GLMGS.Ippt) 

PARAMETER  (KSIZEslOOO) 

REAL  NO.  MU.  GAMLRG.  GAMSML.  N.  ktqKKSIZE  JCSEZE) 

COMMON  /BLOCKl/  NO.  MU.  GAMLRG.  GAMSML.  klQ) 

COMMON  /BLOCK4/  EQ20TI 
C 

C  Compute  the  equation  labelled  20  at  equation  S 
C 

TPPT  »  THUME  +  ITNCR 
EQ20T2  =  EXP(-GAMSML  •  TPRIME) 

EQ20T2P2  =  EXP(-GAMSML  •  T) 

EQ20T2  =  (EQ20T2  -  EQ20T2P2)  /  GAMSML 
EQ20T2N2  =  (EXP(-GAMSML  •  TPPT)  - 

*  EQ20T2P2)  /  GAMSML 
EQ20T3  =  CTPRIME  -  T) 

EQ20T3N2  =  CIPPT-T) 

EQ20T4  =  NO  *  EQ20T2 
EQ20T4N2  =  NO  •  EQ20T2N2 
EQ20T5  =  NO  •  (T  -  TPRIME) 

EQ20T5N2  =  NO  *  CT  -  TWT) 

EQ20p  =  EQ20T1  *  (EQ20T2  +  EQ20T3)  +  EQ20T4  -  EQ20T5 
EQ20N2  =  EQ20TI  •  (EQ20T2N2  +  EQ20T3N2)  + 

*  EQ20T4N2  -  EQ20T5N2 
EQ20  =  (EQ20p  +  EQ20N2)  /  2.0 
IF(((int(t*(1.0/tincr))).GTJCSIZE).OR. 

*  (int(qffime*(l.(ytincr)).gtksize))THEN 
WRITE(6.*)TNCREASE  SIZE  OF  KSIZE  VARIABLE’ 

CALL  EXIT 

END  IF 

ktlp(int(t*(1.0Aincr)).int(tpriine*(l.(Vtincr)))  ^  eq20 
RETURN 
END 
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A  STUDY  OF  THE  DIFFERENCE  IN  PULSE  LINESHAPE 

FROM  SF  AND  ASE 

Amplified  spontaneous  emission  (ASE)  and  superfluorescence  (SF)  can  be 
distinguished  theoretically  (Ref.  C-1)  and  experimentally  (Ref.  C-2).  In  this  section  we 
discuss  the  differences  between  ASE  and  SF  pulse  shapes  calculated  in  the  semiclassical 
approximation.  From  the  work  of  Trammell  and  Hannon  (Ref.  C-3)  we  obtain  the  pulse 
intensity  for  an  initially  Lorenzian  pulse  after  passage  through  an  amplifying  medium  of 
length  /. 

Ip(t)  =  e'^*  Lt  (2e-^‘  - 1 )]  e'^  . 

with 

K  -27tX^fgAn*(t) 

0"  (a+l)(a+l)  ’ 

where 

X  is  the  wavelength  on  resonance 

f  is  the  recoilless  fraction 

g  is  the  coupling  constant  (=  1  for  present) 

An*  is  the  inversion  density  at  t  ==  0 

a  is  the  internal  conversitm  coefficient 
a  is  the  inhomogeneous  broadening  parameter. 

From  the  work  of  Bonifacio  and  Lugiato  (Ref.  C-4)  we  obtain  the  pulse  intensity 
for  an  SF  emission  which  is  given  by 


where 
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ti)  =  1/2  tR  /n  N 

N  is  the  total  number  of  cooperative  nuclei 

^  (l  +  a)(l  +  a) 

XR  =  k  (g2  p)  =  87c 
To  is  the  natural  lifetime 

p  =  N/V  is  the  density  of  radiators  in  the  cooperation  volume. 

In  making  the  comparison  between  Is(t)  and  Ip(t)  we  assume  the  nominal  values  of: 

To  =  I  s, 

p  =  io22  cm-3 

X  =  10~8  cm 

/  =  1  cm 

f  =  10-1 

a  =  10 

a  =0 

V  =  3x10-8  cm3 
T2'''  — >  00  . 

Figure  C-1  shows  the  variations  of  the  ASE  pulse  Ip(t)  as  a  function  of  the 
inversion  An*  for  values  of  An*  from  10^2  cm-3  to  cm"3  and  p,  =  0.  Figure  C-2 
shows  the  variation  in  the  pulse  for  p  =  0, 1, 10, 10^  and  An*  =  1022.  Figure  C-3  shows 
the  SF  pulses  IsfO)  as  a  function  of  N  =  nV  and  Fig.  C-4  shows  the  SF  pulses  as  a 
function  of  p.  To  perform  this  calculation  both  tr  and  Td  were  modified  by  the  factor  0/p 
as  described  in  the  text  to  take  into  account  the  reduction  in  the  effective  cooperative  volume 
for  large  values  of  p. 
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X  it) 


Figure  Ol-  The  ASE  pulse  as  a  function  of  time  for  different  Initial  Inversion 
densities  In  cm-3,  (a)  *  lO**.  (b)  =  10*1 ,  (c)  5  x  1020.  (d)  lO^o.  (e)  5  x  iQi®. 
(f)  lO'**,  (g)  10^®,  and  u  s  0  and  r  =  1  s-’  for  all  cases. 


TIME  (s) 


Figure  C-2.  The  ASE  pulse  as  a  function  of  time  for  different  values  of  u 
(a)  p  s  0,  (b)  p  s  0.1  cm-'',  (c)  p  s  i  cm-'',  (d)  p  =  10  cm-'', 
and  the  Initial  Inversion  An*  s  10^^  cm-3  and  r  s  i  s. 
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METHOD  OF  ANALYSIS  OF 
THERMAL  TRANSFER 

Dl.  THE  FINITE  ELEMENT  SIMULATION  TECHNIQUE 

The  finite  element  simulation  technique  (FEST)  has  been  described  in  Refs.  D-1, 
D-2,  and  D-3  and  has  been  used  to  treat  several  different  types  of  problems.  We  will 
therefore  only  sketch  the  general  procedure  here. 

FEST  applies  the  laws  of  conduction  of  heat  and  mass  actitm  to  a  system  composed 
of  any  number  and  variety  of  materials  in  the  following  way.  The  system  is  divided  into  a 
number  of  convoiient  cells,  each  with  a  definite  density  p,  thermal  conductivity  K,  specific 
heat  Cp,  and  diffusivity  K.  The  temperature  in  each  cell  is  specified  initially  and  then  a  time 
interval  is  chosen  over  which  conduction  is  m  take  place. 

Heat  conduction  during  the  time  interval  At  takes  place  from  one  cell  to  another 
according  to  Fourier’s  law  of  heat  transfer.  For  cells  (i)  and  (i  -  1)  this  is  given  by 

AQ(i  -  1.  i)  =  K  (i.  i  -  1)  cr,  -  T^,)  At 

AX(i,  i-l)  » 

where  A(i,  i  -  1)  is  the  area  across  which  the  transfer  takes  place.  K(i,  i  -  1)  the  effective 
tiiermal  conductivity  between  cells,  which  can  be  calculated  firom  the  thermal  conductivities 
of  the  two  cells,  AX(i,  i  -  1)  the  separation  of  the  cell  center,  Ti  -  Ti  -  1  the  temperature 
difference,  and  AQ(i,  i  - 1)  is  the  heat  transferred  between  cells  i  and  i  -  1  during  the  time 
interval  At 

The  temperature  change  in  the  ith  cell  after  this  transfer  of  heat  between  the  ith  cell 
and  the  two  neighboring  cells  (i  - 1)  and  (i  + 1)  is 

AT.  =  [AQ(i-l.i)-AQ(i.i+l)](CpPjV.)  ,  (D-2) 


D-3 


where  Vi  is  the  volume  of  cell  i.  The  conductance  is  expressed  as  endi  =  Kj  Ai/AXj. 
When  the  medium  changes  at  the  interface  between  i  and  i  +  1,  we  use  instead  the 
expression 

2  endi  cndi+i  /  (endi  +  cndj+i)  . 

This  type  of  program  allows  for  relatively  easy  solutions  to  complicated  multi¬ 
layered  boundary  problems  (Refs.  D-4,  D-5).  Another  medium  can  be  added  without  an 
appreciable  increase  in  complexity  of  the  program,  whereas  in  the  usual  numerical  methods 
of  solving  the  partial  differential  equation  of  heat  conduction  this  is  not  the  case. 

In  addition  to  the  passive  heat  transfer  due  to  thermal  diffusion,  one  can  take 
account  of  heat  sources  in  the  cell  (i)  in  a  natural  way  by  including  in  addition  to 
A(5(i  -  1,  i)  the  term  AQ(i)  =  H(i,  t)  At  which  gives  the  heat  change  due  to  a  reaction 
occurring  in  the  cell. 

D2.  STABILITY  AND  RESOLUTION 

In  any  computational  scheme  a  certain  time  interval  At  must  elapse  during  which 
heat  flows  uninterruptedly  between  cells.  If  the  At  chosen  is  too  large,  then  the  equilibrium 
point  win  be  passed  and  oscillations  in  the  confuted  temperatures  wiU  lead  to  meaningless 
results.  This  is  discussed  in  detail  by  Carslow  and  Jaeger  (Ref.  D-4)  and  Davids  and 
Berger  (Refs.  D-2  and  D-3).  The  well  known  condition  which  insures  against  such  a 
blow-up  is 

XMOD  =  KAt/Ax2  <  1/4  (D-3) 

where  Ax  is  the  effective  linear  dimension  of  the  ceU  in  question  and  k  the  effective  thermal 
diffusivity  between  the  ceUs.  The  above  restriction  on  XMOD  insures  that  during  the  time 
interval.  At,  energy  conservation  is  preserved  and  not  menre  than  the  available  heat  in  a  cell 
is  allowed  to  flow  out 

One  of  the  means  of  reducing  XMOD  below  the  required  value  is  to  increase  the  cell 
size.  This,  however,  tends  to  give  us  less  spatial  resolution  than  may  be  desirable  for  a 
particular  problem.  Furthermore,  it  is  necessary  to  determine  the  number  of  cell  divisions 
that  should  be  made  to  obtain  the  desired  resolution  and  accuracy.  The  general  rule  of 
thumb  is  that  if  S  is  the  dimension  of  the  material,  a  cell  size  of  S/6  gives  a  "reasonable" 
solution.  Each  case  must  be  handled  separately  and  we  have  found  that  sometimes  one  or 
two  cell  divisions  give  sufficiently  good  results. 
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D3.  TWO-DIMENSIONAL  HEAT  CONDUCTION-AXIAL  SYMMETRY 


In  the  heat  conduction  program  with  axial  symmetry,  each  cell  property  which  is 
involved  in  heat  transfer  in  general  is  a  function  of  two  indices,  1  and  J,  corresponding  to 
the  radial  and  axial  directions.  A  representative  cell  is  shown  in  Fig.  D- 1 .  The  geometrical 
cell  properties  are: 


XR(I.J) 

= 

radial  dimension 

ZHa4) 

= 

vertical  dimension 

XS(I4) 

= 

surface  area  in  the  radial  direction 

ZS(I,J) 

= 

surface  area  in  the  vertical  direction. 

Each  cell  is  also  characterized  by  a  thermal  conductivity  K(IJ),  specific  heat  C(I,J),  and 
density  D(I J).  At  each  moment  in  time  the  cell  is  characterized  by  a  temperature  U(I,J), 
and  the  laws  of  heat  conduction  are  used  to  calculate  the  change  in  this  teiiq>erature  over  the 
time  interval  Tl. 


Figure  D-l.  Geometrical  Properties  of  Representative  Cell  (l,J). 


The  various  heat  and  temperature  changes  are  shown  in  Fig.  D-2.  In  the  schematic 
drawing  ZQ(LJ)  is  the  heat  transferred  from  cell  (IJ)  to  cell  (IJ  + 1)  and  ZQ(I,J  -  1)  is  the 
heat  transferred  from  cell  (14  -  1)  to  cell  (14).  The  quantities  XQ(I,J)  and  XQ(I  -  1,J) 
give  the  heat  lost  and  gained  in  the  radial  direction.  The  heat  transfer  quantities  and  the 
resultant  temperature  changes  are  calculated  by  using  the  law  of  heat  transfer  and  mass 
action  as  stated  above. 
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